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The  set  of  programs  and  options  described  in  thisWetSrrt  provide 
a  versatile  new  approach  to  river  basin  planning  for  water  quality  man¬ 
agement.  Altho’.  ,;  :arcr  .quality  relationships  arc  the  primary  considera¬ 
tions  of  the  calculations  performed  here,  the  programs  are  very  general 
and  can  be  manipulated  to  consider  many  water  uses.  Since  water  quality 
strongly  affects  all  uses  of  water,  the  programs  are  designed  to  consider 
most  of  the  legitimate  water  uses.'  The  models  are  capable  of  analyzing 
the  quality-quantity  interrelationships  for  municipal  and  industrial 
water  supply,  flood  control,  fish  and  aquatic  life,  irrigation,  and  rec¬ 
reation.  These  and  other  uses  can  be  considered  either  by  studying  how 
to  manage  existing  situations  or  by  introducing  various  design  alternative- 
concerning  physical  and  structural  components.  As  an  analytical  technkiuc 
for  planning  and  utilizing  river  basins,  this  "River  Basin  Simulation  ‘ 
Program  '  represents  a  major  programming  effort  having  as  its  base  a  sound 
analytical  and  engineering  consideration  of  water  quality.  /  t 


The  programs  offer  a  logical  framework  for  a  sound  technical  analysis 
of  those  aspects  of  a  water  resource  system  that  influence  water  quality. 
As  in  the  case  with  most  useful  and  sophisticated  management  tools,  the 
people  who  commit  themselves  to  using  this  program  will  (1)  find  that 
their  efforts  suggest  useful  alternatives  for  Government  action  and  (2) 
th;y  will  receive  useful  insights  into  the  complexities  of  river  basin 
systems.  Upon  examining  these  and  other  indices  and  distributions,  the 
engineer  on  planner  can  identify  which  alternative  proposal  is  superior. 


The  program  models  the  behavior  of  a  river  basin  with  the  primary 
intention  of  using  Monte  Carlo  Sampling  (simulation)  to  find  the 
occurrence  frequencies  of:  (1)  various  water  quality  indicators  or 

**0j  ^>nt  fluxe?*  ( 2 )  water  supply  deficiencies,  (?)  seasonal  flows, 
and  (4)  reservoir  contents.  The  computations  performed  are  designed 
for  the  users  within  the  Federal  Water  Pollution  Control  Administration: 
hohever,  the  underlying  logic  is  sufficiently  general  to  make  the  program 
useful  to  other  organizations  interested  in  water  resources.  For  example 
irrigation,  flooding,  and  navigation  can  be  studied.  *  ’ 


The  user  prepares  a  sat  of  input  cards  for  the  program  which  serve 
two  functions.  :ne  first  function  describes:  (1)  a  natural  hydrologic 
regime,  (_)  evaporation  and  seepage  conditions,  and  (5)  ambient  water 
quaUty,  The  second  function  specifies  the  plan  which  can  be  co;  prised 

SL!  L!atmqt,aIity10!,j,-CtiVf  (standards) ,  (2)  flow  requirements  and 
diversions,  (J)  waste  loadings  from  tows,  industries,  and  irrigations 
districts,  (1)  various  dans,  and  (5)  (optionally)  a  table  of  the  release 
schedules  to  be  used  at  the  dams. 
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The  program  finds  the  probable  hydrologic  and  water  quality  conse¬ 
quences  of  the  imposition  of  such  a  plan  on  the  unregulated  or  natural 
hydrologic  regime  of  the  basin.  The  j>rogram  permits  the  user  to  inves¬ 
tigate  many  plans  for  the  same  basin. 

Note  that  the  effects  of  various  design  decisions  (or  plans)  are 
measured  in  physical  terms  and  that  there  is  no  facility  in  tinf "program 
to  consider  the  various  costs  and  benefits  involved.  The  economic 
consequences  of  design  choices  must  be  considered  externally  to  the 
program . 

Users  of  the  progra"  can  verify  existing  conditions  and  simulate 
various  management  alternatives.  The  list  of  alternatives  that  can  be 
examined  includes  studying  various  combinations  of  reservoir  sizes, 
reservoir  releases  and  waste  input  schedules.  An  attractive,  major 
feature  is  the  ability  to  study  transfer  of  water  from  another  basin. 

The  study  of  existing  conditions  is  a  means  of  evaluating  the  adequacy 
and  practicality  of  water  quality  standards. 

Answers  are  obtained  in  the  form  of  statistical  information  concerning 
the  probable  consequences  of  the  implementation  of  proposals  under  study. 
Consequences  evaluated  are  water  quality,  streamflow,  and  reservoir  levels 
at  preselected  locations.  The  program  yields  information  including  statis¬ 
tical  moments  and  probability  distributions  of  the  following  quantities: 

1)  Pollutant  concentrations  at  various  prespecified  locations; 

2)  Reservoir  levels; 

3)  Flows  and  truncated  flows  abov'  and  below  prespccificd 
requii-ements. 

For  a  typical  system,  there  is  a  large  amount  of  information  provided 
as  output  from  the  program.  Output  is  labeled  clearly  and  without  ambiguity 
for  the  convenience  of  the  users.  This  wealth  of  information  provides  the 
facts  upon  which  rational  decisions  can  be  based. 

These  programs  are  intended  to  be  an  aid  for  establishing  uniformity 
among  K.'PCA  Project  and  Regional  offices  and  other  agencies,  and  to  make 
statistically  efficient  use  of  hydrologic  or  water  quality  records  that 
exist. 


3 


i 

i 

o 

0 

D 

0 

IS 

n 


n 

Li 

R 


INTRODUCTION 


The  program  offered  here  is  founded  on  well-known  statistical  theory, 
which  is  thoroughly  documented  and  cited  in  the  statistical  and  engineer¬ 
ing  literature.  It  is  presumed  that  at  least  one  member  of  any  team  who 
uses  this  program  is  familiar  with  the  notation,  terminology,  and  under¬ 
lying  mathematical  structure  of  the  theory.  While  intimate  knowledge 
of  the  algebraic  manipulations  is  helpful,  it  is  surely  not  required. 

To  attain  the  necessary  statistical  background,  special  attention 
is  directed  to  the  following  sources: 

1.  Maass,  Arthur  A.,  et  al.,  The  Design  of  Water  Resource  Systems 
(Harvard  University  Press,  Cambridge,  1962),  cf.  chaps.  12^147 


2*  Operations  Research  in  Water  Quality  Management.  Report  of  the 
Harvard  Water  Resources  Group  to  the  United  States  Public  Health 
Service,  1963;  cf.  chaps.  1  and  2. 

3.  Ficring,  Myron  B,  ".Queuing  Theory  and  Simulation  in  Reservoir 
Design,"  Trans  ASCE,  Hydraulics  Division,  127,  1962. 

4.  Fiering,  Myron  B.,  "A  Multivariate  Technique  of  Synthetic 
Hydrology,"  ASCE,  Jour.  Hyd.  Div, ,  90,  September  1964. 

5.  Kendall,  Maurice  G.,  A  Course  in  Multivariate  Analysis,  fllafner 

New  York,  19S7).  “ - < -  ' 


f 
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6.  Fair,  G.  M.  and  Geycv.  J.  C.,  Water  Supply  and  Waste  Water 
Disjoaal^,  (John  Wiley  $  Sons,  TncT,  N.Y.,  1954),  cfTcHps.  4-6. 

7.  Hoel,  P.  G.,  Introduction  to  Mathematical  Statistics.  Ed.  3. 

(John  Wiley  5  Sons,  Inc./T962X. !  ”  '  '  - 

8.  Fruend,  J.  E.,  Statistics,  (Prentice-Hall,  Inc., 

K»  J  * |  1962) « 

9.  Hufschnidt,  M.  M.  and  Ficring,  M.  B,  Simulation  Techniques  for 

^5irj^of_Watcr-ncsourcc  Syst ems ,  PFcsTT ~ 

CamEridgc,  1966X  * 


10. 


Young,, G.  K. 
Residuals," 


and  Pisano,  W.  C»,  "Operational  Hydroloev  Usino 
ASCE,  Jour.  Hyd.  J^TXTW 


t 

i 

! 

i 


) 


4 


Informational  requirements  of  the  program  include  data  on  (1) 
historical  flows  (2)  geometric  location  of  reservoirs,  waste  inputs, 
and  water  users  (3)  background  water  quality  relationships  (4)  evapora¬ 
tion  (5)  magnitudes  of  water  use,  waste  input,  and  reservoir  volumes, 
and  (6)  waste  scheduling  and  reservoir  management  practices. 

Individual  questions  can  be  broad  as,  how  do  present  management 
practices  in  the  basin  influence  water  quality  at  a  state  boundary  or 
how  does  a  change  in  practices  affect  the  quality  at  a  critical  point? 

A  more  specific  question  would  be:  Will  the  flows  in  a  wild  portion 
of  the  river  basin  be  sufficient  to  maintain  a  required  wetted  perimeter 
for  fish  spawning?  Both  types  of  questions  are  important  and  both  can 
be  analyzed. 

Essentially,  there  are  two  elements  in  the  program.  Multisite 
streamflow  simulation  is  designed  to  capture  the  hydrologic  nature  of 
real  river  systems.  This  step  is  designed  to  cause  the  important 
random  characteristics  of  the  historical  records  to  be  duplicated  in 
synthesized  records.  Historical  and  simulated  streamflows  are  routed 
through  the  river  system.  Quality  is  introduced  as  an  operational 
quantity  through  streamflow-water  quality  relationships  as  well  as  by 
defining  the  existing  waste  discharge  practices.  Both  elements  require 
the  full  computing  power  of  the  largest  computers. 

How  big  can  the  river  system  be?  The  river  system  can  have  15 
reaches  where  water  quality  is  critical,  10  confluence  points,  10 
waste  input  points,  and  15  reservoirs.  The  program  can  handle  five 
water  quality  pollutants,  one  of  which  nay  be  degradable,  and  also  can 
consider  dilution  flows  for  water  quality  requirements.  Releases  end 
withdrawals  for  any  purpose  may  be  included  to  study  their  effects  on 
water  quality.  The  unique  feature  of  the  program  is  that  pollutants 
as  well  as  flows  arc  routed  through  the  system.  The  release  schedules 
for  dams  may  follow  current  practice  or  may  be  altered  in  an  attempt 
to  manage  the  system  for  changing  goals  and  needs. 

Each  proposed  design  or  plan  is  statistically  evaluated  in  terms 
of  the  l.ater  quality  and  stror-mflow  characteristics  at  strategic  river 
basin  locations.  A  Mstogr.-.-q  or  frequency  function  of  deficiencies 
from  which  the  mean,  . srliincc,  :nd  skewness  of  the  population  of 
deficiencies  can  bo  v  ;inated,  contains  necessary  design  information 
•bout  the  set  of  deficiencies  which  might  be  obtained  if  the  system 
were  built  or  operated  as  planned. 

For  a  given  river  basin  having  towns,  industries,  rootc&tiennl  areas, 
irrigation  districts,  fish  spawning  reaches,  existing  reservoirs,  and 
perhaps  other  complex  elements,  many  questions  can  be  answered  through 
use  of  the  program.  The  program  relates  the  various  activities  in  the 
basin  to  the  underlying  hydrologic  condition?  which  must  be  satisfied. 


IMPLEMENTATION 


To  use  this  program,  the  following  materials  are  needed.  An 
IBM  360/65  data  processing  system  and  three  classes  of  punched  cards 
describing  option",  controls,  data,  and  design;  these  c  rds  comprise 
the  input.  The  machinery  Is  a  formidable  collection  of  hardware, 
but  it  is  necessary  to  process  the  data  collected  from  a  large  number 
of  sites  in  a  river  basin.  The  capacity  of  the  core  storage  is  nearly 
exhausted  by  the  program.  The  programs  arc  written  in  FORTRAN  IV 
language. 

The  three  classes  of  punched  cards  comprising  the  input  consist 
of  the  following  specifications: 

1)  the  control  cards  which  specify  lengths,  basin  size,  etc., 

2)  data  cards  giving  the  (fixed)  hydrometeorologic  variate  values, 

3)  design  cards  specify  the  set  of  (variable)  design  choices 
to  be  investigated  during  each  run. 

There  are  three  major  options  in  this  package.  One  option  (NMP) 
(Numerical  Multisito  Program)  concerns  the  generation  of  operational 
hydrology.  Another,  (Basin)  is  devoted  to  the  detailed  flow 
routing  necessary  to  model  a  suitably  defined  river  basin.  A  third 
lesser  option  (D),  arbitrarily  called  "DEM**  combines  these  two  concepts 
in  a  one-dimensional  inflow  model;  this  option  is  considered  first. 

1.  Option  D  considers  a  deterministic  or  stochastic  targot  flow 
to  be  met  by  releases  from  a  single  storage  facility  on  a  single  stream. 
Input  consists  of  a  few  control  cards,  followed  by  data  cards  which 
contain  the  hydrologic  record,  the  vectors  of  flow  requirements,  and 
other  miscellaneous  inputs  to  the  model,*  These  are  follov.’ed  by  design 
cards  containing  those  storages  which  are  to  be  evaluated  by  the  program. 

A  sample  single  stream  system  is  displayed  in  Figure  A. 

The  program  deals  in  basic  time  units  or  seasons  of  1,2, 3,4, 6  or 
12  months.  Each  seasonal  flow  requirement  reflects  a  discharge  sufficient 
to  assure  tnat  municipal  and  industrial  water  supply  requirements  arc  met 
or  that  the  maximum  allowable  concentration  of  any  pollutant  is  not 
exceeded  at  a  point  immediately  below  the  dam  unless  all  the  water  in 
upstreau  storage  is  depicted. 

For  eeeh  set  of  seasonal  flow  requirements,  or  dilution  vector, 
the  program  loops  through  all  storage  capacities  provided  in  the  input. 

As  many  as  ten  different  capacities  may  bo  examined. 


•These  miscellaneous  itcas  will  be  davit  with  laicr  in  the  section  on  card 
input  requirements. 
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The  program  allows  for  three  different  nodes  of  reservoir  operating 
policies.  The  first  of  these  modes  always  attempts  to  maintain  the 
downstream  flow  requirement  in  every  season.  The  second  uses  a  set  of 
rule  curves  to  determine  the  outflow  from  the  reservoir  for  each  season. 
This  mode  requires  as  input,  coefficients  for  a  set  of  seasonal  regression 
equations  relating  release  to  storage  and  inflow.  In  addition,  rules 
prescribing  maintenance  of  a  percentage  of  reservoir  capacity  may  be 
Specified. 

J  Stated  in  any  way,  a  rule  can  be  inputted  to  keep  a  reservoir  at 

certain  pool  elevation  during  any  part  of  the  year.  Thirdly ,  the 
'I  program  uses  the  second  mode  with  the  constraint  of  maintaining  down 

stream  flow  requirements. 


_iv. 


Q 
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For  each  combination  of  dilution  vector  and  storage  capacity, 
the  program  routes  the  hydrologic  record  through  the  reservoir  using 
the  prespecified  operating  policy  and  then  tabulates  and  prints 
statistical  characteristics  of  performance.  Finally,  the  program 
generates  univariate  operational  hydrology  (in  accordance  with  the 
Thomas-Ficring  algorithm  which  recursively  produces  synthetic  flows 
in  accordance  with  the  Gaussian  lag-one  Markov  process).  The  operational 
hydrology  is  presumed  to  derive  from  flows  having  a  frequency  function 
that  is  i)  normal,  ii)  square-root  normal,  iii)  log  normal.  The  program 
specifies  the  appropriate  choice  from  the  input  data.  (10). 

In  addition  to  synthetically  generated  (operational)  hydrologic 
inputs,  the  program  can  generate  a  random  component  for  each  season's 
flow  or  dilution  requirement.  The  requirement  in  any  season  is  presumed 
to  vary  normally  about  its  average  value  and  the  input  includes  the 
season  mean  and  the  seasonal  standard  deviation  Aj .  From  these  and 
a  normally  distributed  standardized  sample  deviate,  a  random  number 
u  (0,1)  N  the  program  calculates  Rj,  the  target  flow  for  season  j  from 

Rj  -cj  -u  A,  £0. 

To  run  the  model  using  a  set  of  deterministic  values  for  the  dilution 
vectors,  it  is  necessary  only  to  setAj  -  0 


y 
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If  desired,  seasonal  diversions  may  be  considered  at  the  reservoir  site 
before  any  releases  for  downstream  users  ore  r.ndc .  This  allows  for 
evaporation  losses,  seepage,  or  direct  pumpage  out  of  the  reservoir 
into  another  basin.  Any  addition  to  the  system  below  the  reservoir; 
c.g.,  irrigation  return  flows,  may  he  accommodated.  This  is  accomplished 
by  prefixing  the  seasonal  diversion  with,  a  minus  sign.  As  in  the  case 
of  the  downstream  flow  requirements,  a  random  component  may  be  added 
to  the  seasonal  mean  diversion. 


The  routing  and  tabulation  of  system  performance  under  synthed  '.ally 
generated  variables  (operational  hydrology)  is  effected  identically  Lo 
the  historical  record.  When  all  dilution  vectors  are  exhausted,  the 
program  searches  the  input  tape  for  a  new  site  to  investigate.  Should 
none  be  available,  the  job  terminates. 

This  option  pertains  to  a  single  stream;  the  multi-stream  case  is  studied 
by  the  dual  usage  of  Options  NMP  and  B  as  discussed  below. 

2.  Option  NMP  deals  with  multivariate  distributions  of  hy4?or»«tc<»i’ologic a  1 
data,  with  their  correlation  structure,  and  with  the  generation  of  long  traces 
of  operational  hydrologic  records  using  parameters  extracted  from  the  observed 
record.  The  basic  tine  unit  of  the  generating  algorithm  is  subject  to  input, 
control,  but  is  limited  to  intervals  of  1,2, 3, 4, 6  or  12  months.  See  (10) 

for  operational  details. 

The  input  deck  for  NMP  consists  of  several  control  cards  and  a  large 
number  of  data  cards  giving  the  variate  values. 

3.  Option  B  is  a  generalization  of  the  flow  routing  performed  in 
Option  D.  Instead  of  investigating  the  performance  of  a  system  consisting 
of  a  single  structure  and  a  single  use  point,  this  option  can  consider 
any  geometrical  configuration  of  up  to  15  potential  reservoir  sites,  10 
effluent  discharge  sites,  25  use  or  test  sites,  and  5  pollutants.  Stream- 
flows,  effluent  loads  and  releases  from  reservoirs  are  perfectly  and 
instantaneously  mixed  and  arc  routed  downstream  from  point  to  point  in  the 
system. 


The  entire  routing  of  the  basin,  from  the  headwaters  to  the  mouth, 
is  accomplished  during  each  time  frame.  The  month  is  the  length  of  the 
time  frame  considered  by  this  program.  Thus  only  those  river  basins 
whose  strcamflows  have  an  average  travel  time  of  a  month  or  loss  can  be 
treated. 

This  option  performs  (logically)  much  like  the  simpler  Option  0. 

A  notable  exception  is  that  the  synthetic  variables  ( s t r o am f 1 ows)"  are  not 
generated  by  Option  P  because  of  the  time  and  storage  limitations;  particu¬ 
larly  the  latter.  To  examine  behavior  using  historical  flows  as  inputs, 
the  program  expects  the  observations  to  be  on  cards  or  magnetic  tape;  for 
long  operational  traces  the  program  will  generally  be  run  following  Option  NMP . 
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A .  Specification  of  the  Gcoraet ry  of  the  System 

This  option  operates  on  any  configuration  of  50  or  fewer  points. 

By  definition  each  point  or  coordinate  is  a  location  of  utilitarian 
or  geometric  interest  and  five  types  of  points  are  considered  in  the 
program: 

1.  Test  point  or  use  point;  a  water  supply  demand  or  diversion 
is  treated  as  a  test  point. 

2.  Waste  discharge  point. 

3.  Damsite. 

4.  Junction  point;  the  confluence  of  two  (and  only  two)  streams. 
The  program  considers  all  iunction  points  to  be  test  points. 

5.  End  point. 

A  schematic  diagram  of  the  basin  and  the  interconnections  among  its 
components  is  essential  for  preparation  of  the  card  input.  The 
following  set  of  rules,  and  the  related  example,  demonstrate  the 
order  end  format  of  the  input  deck.  The  sample  basin  is  displayed  in 
Figure  B.  River  miles  to  all  points  arc  shown  in  a  coordinate  system 
approximately  equivalent  to  that  of  the  STORET  System,  but  it  will 
be  seen  that  the  exact  mileages  and  geographic  coordinates  are  far 
less  important  than  the  relative  position  of  each  point  with  respect 
to  its  neighbors.  Development,  in  the  basin  is  indicated  by  the  presence 
of  dams,  municipalities,  and  industries;  the  drainage  area  corresponding 
to  each  reach  of  the  several  streams  is  cross-hatched. 


-,5r.v 
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2.  Continue  numbering  from  the  uppermost  point  on  tho  tributary 
encountered  in  number  1  above.  If  junctions  arc  encountered 
as  the  uppermost  point  is  sought,  select  any  branch  and 
resume  the  numbering  from  its  uppermost  point.  Thus,  point 
4  and  5  arc  located,  but  care  must  be  taken  to  number  each 
junction  only  once  so  that  the  sequence  jumps  over  the  junction 
(point  3)  to  which  a  number  is  already  assigned. 

3.  Continue  downstream,  numbering  each  point,  until  another 
junction  is  encountered.  Thus,  point  6  is  located.  Now 
continue  as  in  number  2  noting  that  the  upstream  junction 
allows  the  choice  of  branch  for  point  7.  The  left  branch 
is  selected,  so  that  points  7  and  8  follow  directly;  point  9 
is  a  junction  whereupon  10  becomes  the  right-hand  branch. 

4.  Points  11  and  12  follow  in  order,  but  a  previous  junction 
(point  6)  now  intervenes  so  that  the  end  point,  downstream 
of  the  junction  at  point  6,  is  assigned  number  13. 

The  following  rules  describe  how  the  program  locates  the  next  point 
to  be  considered,  assuming  that  the  current  coordinate  is  the  ith  from 
among  N  joints. 

1.  A  vector  containing  N  elements  is  punched. 

The  i*"  element  is  the  coordinate  of  the  (i+l)st  point, 
this  necessarily  being  the  next  downstream  point. 

2.  A  junction  point  appears  twice  in  this  vector.  That  is,  it 
is  downstream  of  2  distinct  coordinates. 

3.  Starting  points  in  each  branch  do  not  appear  because  they  arc 
not  downstream  from  any  point. 

4.  The  clement  corresponding  to  the  terminal  point  is  assigned 
zero.  Note  that  this  is  not  necessarily  the  element 
because  it  is  possible  that  the  terminus  will  not  be  the  last 
numbered  or  N*"  coordinate. 

All  junctions  define  an  additional  vector  whose  elements  are  the  coordinates 
of  the  starting  or  uppermost  points  of  the  tributary  encountered  above. 

Referring  to  Figure  C,  tho  vector  cf  next  points  is  (2,3,6,5,3,13,8,9,11,9, 

12,6,0].  ”  ~  “  ~ 

The  undorscc:  !  elements  arc  junctions,  and  appear  twice.  Starting  points 
(1,4,7,10]  do  not  appear,  and  tho  Nt5'  or  13th  value  is  zero.  The  vector 
of  starting  p  hits  for  each  junction  is  (4,7,10], 

Starting  point  (1),  being  the  initial  or  leginning  point  in  the  system,  is 
defined  by  other  input. 
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'1  he  program  considers  the  following  transformed  system  of  streams 


Flows  and  effluent  flux  ivies  arc  referenced  to  a  set  of  discrete 


points  which  represent  continuous  tones  or  reaches  of  the  prototype.  Mule 
this  is  a  departure  from  reality,  the  errors  introduced  by  concentrating 
the  processes  of  a  basin  into  a  set  cf  point  processes  are  almost  surely 
within  the  uncertainties  and  errors  propagated  by  inaccuracy  in  the  data. 


k 


Rules  for  establishing  a  coordinate  system  for  system  specification 

arc: 


1.  Start  from  the  uppermost  point  cn  any  branch  and  number  the 
points  consecutively  down  to  and.  including  the  first  junction 
point.  Consider  the  western-  >st  branch  of  Figure  C  to  be  the 
stalling  branch,  so  that  points,  1,  2,  and  3  are  located  directly; 
th?  process  is  interrupted  at  point  3  because  it  is  a  junction 
point . 
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B.  Relation  of  coord) nates  to  flow  gages. 

As  many  as  50  gaging  stations  may  be  used;  Option  KMP  is  limited 
to  25  sites  so  that  if  long  operation-1  traces  arc  desired  Option  I? 
is  constrained  to  25  gages.  A  vector  of  N  elements  is  punched,  the 
1th  element  being  the  gage  number  to  \  hi eh  the  ita  point  or  coordinate 
is  referenced.  The  gaging  stations  are  not  included  in  the  coordinate 
system.  An  additional  vector  provides  drainage  area  factors  so  that  the 
flow  at  any  gage  can  be  scaled  linearly. 

Two  options  arc  available  to  scale  the  gaging  station  flows.  For 
the  first  option,  the  i*^  drainage  area  factor  multiplied  times  the  flow 
of  the  referenced  gage  defines  the  flow  between  the  i^1  point  or  coordinate 
and  the  (i - 1 )  coordinate  or  the  point  directly  upstream.  This  option  is 
used  when  the  gaging  stations  are  located  on  the  tributaries.  If  the 
second  option  is  specified,  then  the  i^  drainage  area  factor  multiplied 
times  the  gage  flow  defines  the  total  inflow  into  the  coordinate. 

c  Testing  Different  Plans 

Because  the  input  dec!:  is  large,  it  is  useful  to  be  able  to  evaluate 
additional  plans  using  the  same  system  configuration  by  c'.  :g  th> 
original  design  cards.  Under  input  control  they  may  be  entirely  re¬ 
placed,  slightly  altered  (a  few  cards  changed),  or  left  intact.  Each 
additional  plan  comprises  a  run.  Each  run  may  use  new  hydrologic  inputs 
or  may  be  instructed  to  leiocate  those  same  flow  data  used  in  earlier 
runs  so  that  system  performance  can  be  compared  under  the  same  stochastic 
regime. 


D.  Wastes  or  Pollutants. 


Five  waste  variables  may  be  considered  simultaneously  in  the  p.ogn 
They  may  be  introduced  at  up  to  10  points.  The  pollutants  are  characterized 
by  seasonal  means  and  standard  deviations.  A  stochastic  component  is  added 
to  the  expectations  of  both  waste  flow  and  waste  weight,  the  magnitudes 
being  proportional  to  the  specified  respective  standard  deviations. 

Three  different  types  of  pollutants  may  be  considered.  The  first 
arc  conservative  pollutants  having  no  counterparts  in  the  natural  or 
ambient  system;  for  example,  industrial  wastes  such  ns  chromates.  These 
wastes  arc  imposed  upon  a  system  having  ?  ro  natural  concentrations. 

The  second  arc  conservative  pollutant:.;  having  counterparts  in  the  natural 
system.  Typical  examples  would  be  api cultural  and  industrial  wastes 
that  include  pollutants  such  as  total  dissolved  solids,  boron,  chlorides, 
and  sulfates,  l.astly,  n  degradable  pollutant  such  as  t.0.1).  nay  be 
considered.  For  this  type  i  r  waste  the  background  or  natural  concentra¬ 
tion  nay  be  included. 
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The  program  routes  up  to  5  wastes  by  attempting  to  maintain  target 
(maximal)  concentrations  at  the  several  test  and  use  points.  In 
addition,  an  unspecified  pollutant  or  set  of  pollutants  can  impose  a 
minimum  flow  requirement  (determined  externally)  which  the  system 
seeks  to  maintain.  Dissolved  oxygen  standards  are  commonly  treated 
in  this  fashion. 


Only  one  degradable  waste  may  be  specified,  and  if  so,  it  is  defined 
to  be  waste  number  one.  For  the:  other  wastes  (type  3  or  type  2)  the 
order  is  not  important.  However,  once  an  order  is  established,  it  must 
be  the  same  for  all  waste  sites.  Furthermore,  the  same  order  must  be 
maintained  for  the  water  quality  standards  or  targets  at  the  test  sites. 

For  each  conservative  pollutant  (type  2),  the  background  concentration 
of  the  unregulated  inflow  at  each  coordinate  in  the  system  is  calculated 
by  one  of  two  options.  The  first  option  uses  log  concc-ntral i on  versus 
log  flow  regression  equations,  'lhe  parameters  of  these  regression 
equations  are  required  as  input.  They  are  determined  fro:,  the  analysis 
of  the  historical  strcemflow  and  cjuaJJjty  observations.  Also  required 
as  input  are  the  uppcr'ancl  lower  concentration  limits  for  each  relation¬ 
ship.  These  limits  are  determined  from  observation  of  the  historical 
data. 


For  the  second  option,  streanflow  and  water  quality  monitoring 
stations  arc  simulated  simultaneously  in  option  In  option  l’>  the 

concentration  at  each  site  is  computed  from  the  simulated  water  quality 
at  the  monitoring  stations  by  linear  regression.  The  parameters  of  these 
relationships  arc  required  as  input.  These  arc  determined  from  the 
analysis  of  historical  quality  observations. 

If  a  degradable  waste  is  specified,  the  following  information  is 
required: 

1) .  Natural  cr  background  concentration  for  each  point;  this  is 

the  load  introduced  from  the  unregulated  flow. 

2)  Mean  seasonal  temperature  (°C.). 

5)  Fiver  miles  from  each  coordinate  to  next  coordinate. 

4)  The  progra  :  computes  a  velocity  V  in  miles  per  day  from  the 
relationship  V  *  a  Q* ,  where  p  and  b  are  required  ns  inputs 
for  each  site, 

5)  I'ooxy,  .  nation  coefficients  (?0*C.)  for  each  point  ami  for  each 
dam  site. 

In  Figure  C,  the  waste  points  are  numbered  (1,2, >,-!). 


E.  Dam  Sites 


IS 


For  each  of  15  or  fewer  sites,  the  capacity,  minimum  pool  storage, 
and  starting  active  storage  arc  required.  Perfect  mixing  is  assured 
in  reservoirs  and  is  accomplished  before  releases  are  made.  As  an 
option  the  combined  reservoir  losses  resulting  from  evaporation  and 
seepage  can  be  considered. 

In  Figure  C,  the  dams  arc  numbered  (1,2,3, 4)  os  shown. 

F.  Reservoir  withdrawal  and  release  options 

A  series  of  choices  may  he  selected  to  operate  any  reservoir  in 
the  system.  Any  reservoir  may  operate  by  combination  of  these  choices; 
however,  the  priority  of  the  options  is  fixed  for  all  reservoirs  by  the 
program.  During  each  time  frame  the  program  performs  ar  initial  routing 
using  only  choice  1  and  choice  2.  If  supplemental  flow  needs  arise 
in  the  basin  and  if  either  of  the  latter  two  choices  have  been  specified, 
then  the  program  performs  several  additional  routings  using  all  the 


specified  choices  to  satisfy  flow  requirement  needs. 


1.  Explicit  (pipeline)  withdrawal.  Seasonal  quantities  of  water 
and  its  quality  may  ’  e  withdrawn  (pipeline)  from  any  reservoir  and  mixed 
and  blended  at  a  test  site.  The  addition  of  water  is  jnr.de  regardless 
of  the  state  of  the  receiving  test  site.  Consider  the  following: 


Shown  in  both  applicable  examples  is  n  reservoir  and  a  test  site. 

Rule:  The  streanflow  at  the  test  site  rist  not  enter  the  dan  site. 

Yhc  dam  must  be  encountered  prior  to  the  test  site  during  the  routing. 

The  typical  operation  for  a  rerrvoJr  using  this  choice  is  the  following: 

1)  Yhc  current  inflow  and  start  of  the  season  storage  arc  added 
and  the  quality  of  the  contents  blended. 

2)  Yhc  evaporation  loss  or  rainfall  addition  is  then  considered 
and  the  quality  te- con«t Itutcd. 
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3)  The  required  diversion  is  then  tdthdrawn  from  the  total  assets. 

4)  If  no  other  rclca.se  option  has  been  specified,  then  continuity 

is  cheeked. 

6)  When  the  test  site  is  encountered  in  the  routing,  the  diversion 

is  added  to  the  flow  at  the  test  site  and  the  qualities  blended. 

The  seasonal  quantities  to  bo  withdrawn  are  required  as  input  for  each 
reservoir.  Required  also  is  a  vector  relating  the  locations  of  the 
receiving  test  sites  to  that  of  the  reservoirs. 

2.  Operating  rules.  Any  reservoir  may  be  operated  by  the  following 
seasonal  linear  regression  relationships: 

°ifcAi+Bi*Vc!si  *  =  seoson 

where  Dj,Qj,  end  Sj  are  the  seasonal  draft,  inflow,  and  storage,  respectively. 
Aj.Bj,  and  Cj  ere  ii  anted  for  each  reservoir. 

A  reservoir  nay  be  maintained  at  a  certain  pool  level  fcy  setting 
fcj  and  C.  equal  to  j'ero  and  letting  A.  equal  the  percentage  of  total 
capacity1implicd  by  that  pool  elevation.  Aj  is  input ed  as  a  negative 
decimal.  Recreational  and  flood  control  pool  elevation  rules  ray  bo 
expressed  in  this  manner. 

To  nap  the  operation  of  existing  reservoirs  into  a  system,  perform 
regression  analysis  on  the  observed  inflows,  storage  volumes,  and 
releases.  The  computed  coefficients,  Aj ,  B  j ,  and  Cj  may  then  be 
input cd  to  the  program.  Nolo  that  the  effects  of  evaporation  for  the 
existing  reservoirs  is  already  implicit  in  these  coefficients. 

Minimum  flow  release  for  fish  protection  or  navigation  may  be 
accomplished  by  setting  In  and  Cj  equal  to  aero  and  letting  Aj  equal 
the  minimum  required  release. 

Typical  operation  is  t lie  following: 

1)  The  current  inflow  and  start  of  tliv  season  storage  arc  added 
and  the  quality  of  the  co. -.tents  Mended. 

2)  The  evaporation  loss  or  rainfall  addition  is  then  considered 
•nd  the  quality  re-constituted. 

3)  If  an  explicit  diversion  is  specified,  it  is  withdrawn. 

4)  Release  computed  fro  >  the  seasonal  rule  is  then  drafted. 

5)  Continuity  is  checked. 


3.  Conditional  pipeline  diversions.  Any  test  site  nay  receive 
additional  waters  to  satfsfy  a  minimum  flow  requirement  and/or  water 
quality  standards.  This  choice  differs  from  the  first  in  that  the 
conditions  at  the  test  sites  arc  first  investigated  to  sec  if  supple¬ 
mental  waters  are  needed.  Consider  the  following  example: 


Show  in  the  example  is  a  test  site  and  a  reservoir.  The  test  site 
represents  a  city  having  a  water  supply  flow  requirement.  Th  :  city 
obtains  its  water  supply  from  the  runoff  above  it  and  draws  upon  the 
reservoir  only  when  the  natural  source  is  deficient. 

Rule:  The  dam  site  must  be  encountered  prior  to  the  test  site 
during  the  routing. 

Typical  operation  is  the  following: 

1)  Explicit  withdrawal  is  considered  first. 

2)  Release  defined  by  the  operating  rule  (choice  2)  is  then  drafted. 
In  absence  of  choice  1  and  choice  2,  the  reservoir  contents  ere 
checked  for  continuity.  Excesses  are  released  and  proceed  down¬ 
stream  of  the  reservoir. 

J)  The  supplemental  requirement  is  then  withdrawn  from  the  reservoir. 

4)  h’hen  the  test  site  is  encountered  in  the  routing,  the  diversion 
is  added  to  tire  flow  at  the  test  site  and  the  qualities  arc 
blended. 

A  vector  relating  the  locations  of  the  test  sites  receiving  supplemental 
flows  to  that  of  the  reservoirs  is  required  as  input. 

4.  Conditional  channel  re I  cusps .  Any  test  site  downstream  from  a 
group  of  reservoirs  f  <lrt)  can  draw  from  the  storage  of  these  reservoirs 
when  a  need  arises  to  satisfy  a  mini: .urn  flow  reqv * •  cr.ant  and/or  water 
quality  standards.  The  releases  are  iwule  directly  into  the  channel  and 
flow  down  to  the  reed  point.  Specifying  second  tier  reservoirs,  that  is, 
reservoirs  in  series,  is  not  allowed.  Hows  released  to  satisfy  headwater 
needs  are  considered  in  computing  the  supplemental  requirements  of  the 
downstream  test  sites. 
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Three  alternatives  are  available  to  distribute  a  need  among  a  scries 
of  reservoirs:  1)  Demand  is  evenly  split,  2)  Demand  is  proportionally 
split  by  the  r.vrilnbJc  storage  in  each  of  the  reservoirs,  3)  Demand 
is  satisfied  by  sequentially  depicting  a  series  of  reservoirs.  The 
order  of  depletion  is  specified  by  input. 

Consider  the  following: 

T1  & 


T3 


Shown  in  the  example  arc  three  reservoirs  and  three  test  sites  having 
minimum  flow  requirements,  Reservoirs  two  and  three  have  recreational 
and  flood  control  pool  elevation  rules.  Test  site  one  may  receive  a 
conditional  pipeline  diversion  from  reservoir  one.  Test  sites  tv;o  and 
three  nay  draw  upon  reservoirs  two  and  three  respectively  when  the  flow 
requirements  are  not  satisfied. 

The  flow  requirement  at  test  site  two  is  compared  with  the  local 
inflows  into  that  site  plus  any  releases  wade  to  maintain  the  pool  level 
of  reservoir  two.  If  the  flow  requirement  is  not  satisfied  then  additional 
waters  arc  released  from  reservoir  two. 

For  test  site  three,  the  flew  requirement  is  compared  with  the  sum  of: 

<0  all  local  flows  between  the  dens  and  the  test  site;  M  the  conditional 
diversion  withdrawn  for  test  site  one;  c)  ihc  sun  of  all  releases  made  to 
cither  maintain  continuity  or  satisfy  the  pool  level  rules;  and  d)  any 
additional  release  required  to  satisfy  test  site  two.  If  a  deficiency 
occurs,  then  this  amount  is  released  by  reservoir  three. 

5.  Combination  of  choices  and  priority.  Any  reservoir  may  be  operated 
by  any  comFinat ion  of  the  st;  too  ~ch.  ices;  "how *er,.  the  priority  of  releases 
is  fixed  for  all  reservoirs  by  the  program  The  following  is  the  priority 
followed:  1)  explicit  pipeline  diversions,  2)  operating  rule  release- con¬ 
tinuity,  3)  conditional  pipeline  diversions,  and  d)  conditional  channel 
releases, 

f>.  Vast  i*c*»v»s.  All  test  points,  including  all  junction  points, 

vvqXi  "  |  sx-icnts  of  a  flow  target  and  seasonal  quality  targets  for 

esen  *  t c  variable.  Three  different  type*  of  flow  targets  may  be  considered. 
The  fi  are  instv*  ?  .  4n»r.o.a  fitw  icquirvents.  Co.  .on  examples  are 
t-in  i  *  Vet  :w.  ini  lining  dissolved  oxygen  standards,  for  ensuring  fish 

pas  .  ;  f-pe-niug ,  ?nd  fur  enhancing,  recreational  areas.  The  second  are 


n 

!  j 


f 

■! 


minimum  flow  withdrawals;,  that  is,  ..mounts  of  water  that  arc  exhausted 
from  the  system.  Often  cities  receive  their  water  supply  by  diverting 
waters  from  another  basin.  Da.'i.v-le ucics  for  both  arc  defined  as  missing 
or  not  maintaining  the  flow  target.  The  third  type  is  an  instream  flow 
requirement  but  differs  from  the  first  in  that  the  deficiencies  are 
defined  as  excesses  above  the  target  flow.  Here  the  flow  requirement 
may  he  the  channel  capacity  of  a  reach  where  flood  damages  can  be  roughly 
assessed.  The  program  does  not  attempt  to  minimize  these  excesses;  it 
only  tabulates  them.  Figure  C  shows  four  test  points  in  the  system. 
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CARD  INPUT  1-OIi  OPTIONS  !>,  K.TP,  AND  D 


Input  to  the  options  is  arranged  according  to  the  following  descriptions, 
where: 


The  symbol  * 


denotes  that  a  series  of  cards  as  described  may 
be  needed. 


The  notation  P 


indicates  "right  hand  justified", 
quantity  so  described  must  appear 
possible  tc  the  right  of  its  data 


i . e . ,  any 
as  far  as 
field. 


The  symbol  . 


indicates  that  a  decimal  point  must  appear  in 
the  field,  and 


The  combination  P 


indicates  that  the  value  is  right-hand  justified 
integral,  but  may  have,  if  necessary,  a  decimal 
point  punched  in  any  column  to  override  the  pro¬ 
grammed  decimal  which  exists  at  the  extreme 
right-hand  boundary  of  the  field. 


The  symbols 


indicate  the  start  of  a  new  card  or*  the  continua¬ 
tion  of  fields  on  the-  same  card. 
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INPUT  FOR  OPTION  D 


CARD 


COLUMNS 


1 

2 


1-72 

1-SR 

6-10R 

11-15R 

16-20R 


*3  1-6R. 

7-12R. 


DESCRIPTION 


Title  of  Study. 

Years  of  historical  record.  (£100). 

Seasons  per  year  (1,2, 3, 4, 6,  or  12). 

Year  in  which  historical  record  begins. 

Calendar  season  corresponding  to  season  1 
of  input.  For  example,  the  USGS  water 
year  begins  in  October,  so  that  if  the 
historical  data  were  obtained  from  Geolog¬ 
ical  Survey  compilation  reports,  this 
parameter  would  be  10. 

Historical  flow  at  gaging  station,  year  1, 
season  1. 

Historical  flow  at  gaging  station,  year  1, 
season  2. 


67-72R.  Historical  flow  at  gaging  station,  year  1, 

season  12. 


1-6R. 


Historical  flow  at  gaging  station,  year  2, 
season  1. 


.  historical  flow  at  gaging  station,  last 

year,  last  season. 

4  1-6R.  Factor  by  which  flows  for  season  one  are 

*  multiplied  to  convert  them  from  gage  to 

dam  site. 


~  W.‘ 


7-12R. 


Factor  for  season  2. 


1-6R. 

7-12R. 


1-5R 

6-10R 

11-15R 

16-20R 

25 

30 

35 

2 


Factor  for  last  season 

Factor  by  which  inflows  at  dam  site  are 
multiplied  to  compute  local  (unregulated) 
flows  between  dam  site  and  use  point  for 
season  1. 

Factor  for  season  2. 


Factor  for  last  season. 

Number  of  seasonal  flow  requirement  vectors. 

Years  of  operational  (synthetic)  hydrology 
to  be  routed  through  the  proposed  reservoir. 

Number  of  reservoir  sizes,  (£10). 

Random  number  generator  trigger. 

Any  odd  integer  greater  than  zero. 

Reservoir  operating  policy  modes, 

1- Meet  downstream  dilution  requirement. 

2- Rule  curves. 

3- Rule  curves  with  constraint  to  meet 
downstream  flow  requirement. 

If  each  synthetic  trace  is  to  be  identical 
for  each  routing,  punch  1. 

If  the  historical  record  routing  is  to  be 
deleted,  punch  1. 

Print  suppress.  Punch  l  wi  desired, 
historical  record  analysis. 

Statistics  of  one-season  misses,  punch  1. 
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Statistics  of  annual  misses,  punch  1. 


Statistics  of  consecutive  seasonal 
misses  (36  seasons  -  maximum  duration) , 
punch  1. 

Seasonal  and  annual  failure  histograms  in 
tabular  form,  punch  1. 

Seasonal  and  annual  failure  histograms  in 
both  tabular  and  graphic  form,  punch  2. 

Seasonal  histograms  of  storage  levels  ir 
tabular  form,  punch  I. 

Seasonal  histograms  of  storage  levels  in 
both  tabular  and  graphic  form,  punch  2. 

Seasonal  statistics  of  flows  past  test 
point,  punch  1. 

Synthetic  trace  analysis 

Statistics  of  one-season  misses,  punch  1 

Statistics  of  annual  misses,  punch  1. 

Statistics  of  consecutive  seasonal  misses, 
(36  seasons  —  maximum  duration) ,  punch  1 . 

Seasonal  and  annual  failure  histograms  in 
tabular  form,  punch  1. 

Seasonal  and  annual  failure  histograms  in 
both  tabular  and  graphic  form  punch  2. 

Seasonal  histograms  of  storage  levels  in 
tabular  form,  punch  1. 

Seasonal  histograms  of  storage  levels  in 
both  tabular  and  graphic  form,  punch  2. 

Seasonal  statistics  of  flows  past  test 
point,  punch  1. 

Yearly,  grand^and  seasonal  statistics  of 
up  to  first  100  years  of  synthetic  trace. 


1-6R. 

7-12R. 

1-6R. 

7-12R. 


55-60R. 


1-6R. 

7-12R. 


1-6R. 

7-12R. 


1-6R. 


Size  of  first  trial  reservoir,  (units 
must  be  consistent  with  streamflow) . 

Size  of  second  trial  reservoir. 


Starting  storage  of  first  trial  reservoir. 
Starting  storage  of  second  trial  reservoir. 


Starting  storage  of  last  trial  reservoir. 


If  reservoir  release  mods  1  is  exercised, 
skip  *10. 

Res ervoir  release  m le  input. 

1)  .  Release  (season)  =  A(scason)*  B(season)* 

Inflow (season)  +  C (season) *Storage (season) . 

2) .  For  a  pool  level  rule^ B (season)  and 

C(season)  must  be  zero.  A(season)  = 
-Percentage  of  capacity,  (decimal) 

3)  .  Rules  are  the  same  for  each  trial 

reservoir. 

Coefficient  A,  season  1. 

Coefficient  A,  season  2. 


Coefficient  A,  last  season. 


Coefficient  B.,  season  1. 
Coefficient  B,  season  2. 


Coefficient  P.,  last  season. 


7-12R. 


Coefficient  C.,  Season  1. 
Coefficient  C,  season  2. 


\  ' 
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11 

1-72 

12 

1-6R 

7-12R. 


Coefficient  C,  last  season. 

Description  of  dilution  requirement  vector. 
Mean  dilution  requirement,  season  1. 

Mean  dilution  requirement,  season  2. 


67-72R. 

Mean  dilution  requirement,  last  season. 

13 

1-6R. 

Standard  deviation  dilution  requirement, 
season  1. 

7-12R. 

Standard  deviation  dilution  requirement, 
season  2. 

67-72R. 


Standard  deviation  dilution  requirement, 
last  season. 


14  1-6R.  Diversion,  season  1.  (evaporation  or  any 

loss  of  storage  at  reservoir  site) .  An 
addition  to  the  contents  of  the  reservoir 
is  implied,  if  the  field  is  prefixed  with 
a  negative  sign. 

7-12R.  Mean  diversion,  season  2. 


IS  1-6R. 


7-12R. 


Mean  diversion,  last  season. 

Standard  deviation  diversion,  season  1. 
Standard  deviation  diversion,  season  2. 


67-72R.  Standard  deviation  diversion,  season  12. 

REPEAT  CARDS  11-15  FOR  EACH  DILUTION 
48  VECTOR. 

RETURN  TO  CARD  6  FOR  A  NEW  SET  OF  RESERVOIRS 
AND  DILUTION  VECTORS. 

RETURN  TO  CARD  1  FOR  NEW  STUDY. 


I 
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INPUT  FOR  OPTION  NMP 


CARD 

COLUMNS 

DESCRIPTION 

1 

1-72 

Name  of  River  Basin 

2 

1-SR 

Years  of  historical  record  (£100) 

6-10R 

Seasons  per  year  (1,  2,  3,  4,  6,  or  12) 

11-15R 

Number  of  gaging  sites  (£.25) 

16-20R 

Years  of  operational  hydrology  or  quality 

21-25R 

Random  Number  Trigger 

Any  odd  integer  excluding  zero 

3* 

1-72 

Name  of  gaging  site  1 

•  •  • 

Name  of  last  gaging  site 

4* 

Flow  or  Quality  Data 

Site  1,  Year  1 

1-6R. 

Quality  or  flow  for  season  1 

7-12R« 

Quality  or  flow  for  season  2 

•  •  • 

•  •  • 

67-72R. 

Quality  or  flow  for  last  season 

«  t  • 

Site  1,  Year  2 

1-6R. 

Quality  or  flow  for  season  1 

•  •  • 

•  •  • 

«iii 

Quality  or  flow  for  last  season 

Site  i,  I'.st  Year 

1-6R. 

Quality  or  flow  for  season  1 

7-12R. 

Quality  or  flow  for  season  2 

•  •  • 


COLUMNS 


CARD 


DESCRIPTION 
Last  site,  last  year. 

1-6R.  Quality  or  flow  for  season  1 


*  *  4 


Quality  or  flow  for  last  season 
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INPUT  FOR  OPTION  D 

* 

Geometry  Cards 
Card  Columns 

1  1-72 

2  4-5R 
9-10R 

14-15R 

19-20R 

24-25R 

*  3 


5 

10 


50 
•  «  ' 

S 


•Supply  card  with  ?  5 


(Sec  following  page  for  typical  input  avrangsicnt) . 

Title  or  identification 

Number  of  coordinates  $50 

Starting  point  coordinate 

Seasons  per  year  (1,2,3,4,6,12) 

Number  of  gaging  stations  <50,  but  f25  if  flows 
arc  generated  previously  in  Model  .NMP . 

Note:  If  both  stream  flows  and  water  quality 
(conservative)  have  been  generated,  punch  the 
sum,  <25. 

Number  of  complete  runs  (number  of  plans  to  be 
investigated) 

Up  to  five  cards  may  be  needed 
Type  of  point  or  coordinate 

1  -  Test  point  sl5. 

2  -  Waste  point  <10 

3  -  Dam  site  <15 

4  -  Junction  point  <10 

5  -  End  point  fonly  one) 

Type  for  coordinate  1. 

Type  for  coordinate  2 

•  «  i 

Type  for  coordinate  10 
Type  for  coordinate  11 

•  i  •  i 

Type  for  last  (Nth  coordinate) 
punch  in  column  1 . 


Geometry  cards  include:  *l-*7 
Control  cards  include  :  8, *9, *10 

Design  cards  include  :  *ll-*24,  47*  -  *49 

Data  cards  include  :  *25-*  46 


OPTION  B 

TYPICAL  INrUT  CARD 
ARRANGEMENT 


Up  to  five  cards  may  be  needed. 
Next  Point  cards 


*4 

4-5R  Coordinate  number  directly  downstream  of  the 

first  coordinate. 

9-10R  Coordinate  number  directly  downstream  of  the 

second  coordinate. 

•  •••  •••• 

49-50R  Coordinate  number  directly  downstream  of  the 

tenth  coordinate. 

*  i  *  * 

4-5R  Coordinate  number  directly  downstream  of  the 

eleventh  coordinate. 


*5 


Last  coordinate  -  punch  zero. 

Up  to  five  cards  may  be  needed. 


4-5R 


Caging  station  to  which  first  coordinate  is 
referred. 


49-50R  Caging  station  to  which  tenth  coordinate  is 

referred. . 

•  •  •  • 

4-5R  Gaging  station  to  which  eleventh  coordinate  Is 

referred. 


•  •••  Caging  station  to  which  last  coordinate  is 

referred. 

A  Up  to  ten  enrds  way  be  needed. 

See  $-30, 

< 

1-10R*  Drainage  .  «  ratio  referring  coordinate  1 

to  jts  i effective  gaging  station.  (This 
ratio  t^v^unlv  eniy  for  the  inflow  between 
coordi  *c  (!)  n.J  coordinate  (1  -  1)3 » 

At  junction  point  the  ratio  accounts  for  Inflow 
below  both  tributary  coordinates. 


t  •  «  •  » 
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41-50R.  Drainage  area  ratio  referring  coordinate  5 

to  its  respective  gaging  station.  (This 
ratio  accounts  only  for  the  inflow  between 
coordinate  (i)  and  coordinate  (i-1)). 


1-10R.  Drainage  area  ratio  referring  coordinate  6 

to  its  respective  gaging  station.  (This 
ratio  accounts  only  for  the  inflow  between 
coordinate  (i)  and  coordinate  (i  -  1)). 


Drainage  area  ratio  referring  last  coordinate 
(This  ratio  accounts  only  for  the  inflow 
between  coordinate,  (i)  and  coordinate 

(i-1)). 


1-10R.  Drainage  area  ratio  referring  coordinate  1  to 

i^s  respective  gaging  station.  (This  ratio 
accounts  for  the  total  inflow  into  this 
coordinate.) 

11-20R.  Drainage  area  ratio  referring  coordinate  2 

to  its  respective  gaging  station.  (This 
ratio  accounts  for  the  total  inflow  into 
this  coordinate.) 


1-10R.  Drainage  area  ratio  referring  coordinate  6 

to  its  respective  gaging  station.  (This 
ratio  accounts  for  the  total  inflow  into 
this  coordinate.) 


.  Drainage  area  ratio  referring  last  coordinate 

to  its  respective  gaging  station.  (This 
ratio  accounts  for  the  total  inflow  into  this 
coordinate.) 

4-5R  Starting  point  of  complimentary  branch  for 

junction  point  1. 


9-10R 


4 


Starting  point  of  complimentary  branch  for 
junction  point  2. 


Starting  point  of  complimentary  branch 
for  last  junction  point  10. 

Note:  Path  of  r.ovcaent  is  conditional 
upon  starting  point. 


Control  Cards 

Card  Columns 


CARDS  BEYOND  THIS  POINT  ARE  NEEDED  FOR  EACH  RUN 


S  1-5R 

10 


IS 


17-20R 


2S 


< 


Number  of  years  of  historical  or 
generated  record. 

Number  of  waste  variables  <5,  (1,2, 3, 4, 
or  5) 

Iv’aste  sites  with  zero  waste  variables 
are  permitted.  ~ 

Logical  tape  unit  for  input  flows.  If 
card  input,  use  5;  if  binary  tape  input, 
use  4.  Card  input  implies  usage  of 
historical  flows.  Binary  tape  input 
tseans  that  a  tape  of  operational  hydrology 
has  been  previously  prepared  in  Option 
KMP. 

Sequential  year  on  the  flow  tane.  reckoned 
from  the  start,  at  which  routing  will 
begin.  The  flow  tape  (binary)  will  be 
rewound  at  the  end  of  each  plan. 

Multi-run  key.  If  single  run,  punch  zero,  if 
multi- run,  punch  zero  for  first  run  and 
one  for  subsequent  runs.  It  is  desirable 
after  the  initial  run  tc  change  only  a 
few  design  parameters  and  proceed  without 
having  to  replace  the  entire  deck  starting 
with  this  card.  If  this  choice  is 
exercised,  include  card  9,  then  skip  to  card 
*48.  If  historical  flows  are  used,  they 
must  be  included  for  each  run  and  after 
the  design  update  cards.  If  this  choice 
is  not  exercised,  then  replace  deck  from 
this  card. 
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30  If  gaging  station  scaling  factors  (*6)  have 

been  prepared  by  choice  A,  punch  zero.  If 
prepared  by  choice  h,  punch  1. 

35  If  any  reservoir(s)  are  to  make  explicit  pipeline 

diversions  to  any  *est  sites,  punch  1. 

40  If  any  reservoir(s)  are  to  make  conditional  pipe¬ 

line  diversions  to  any  test  sites,  punch  1. 

41-45R  Random  number  generator  trigger.  Any  odd 

integer  greater  than  0. 

50  Allocation  alternatives  for  reservoirs  serving 

prc-specified  downstream  test  sites.  Only 
aoplicablc  when  more  than  one  reservoir 
?  .s  been  specified  to  serve  a  downstream 
test  site. 

Punch : 

1- even  split 

2- apportion  by  available  storagefspacc  rule) 

3- sequential  depletion.  Priority  is  speci¬ 
fied  by  entry  order  in  *22. 

55  Mode  of  calculating  conservative  pollutant 

concentrations . 

A.  If  no  conservative  pollutants  are  considered 
in  analysis,  punch  zero. 

B.  If  background  salt  concentrations  arc  to  be 
calculated  by  exponential  regression  equations 
using  streamflow,  punch  1. 

B  •  * 

Cjj  =  Aij^i  S  i-coordinate  index 
where  A- . ,B^ • ,L- . ,Tj •  are  input,  j-pollutant  index 

C.  If  salt  concentrations  at  water  quality  moni¬ 
toring  watering  stations  are  simulated  together 
with  streamflov.s  in  Model  NMP,  then  punch  2. 

Only  one  pollutant  is  permitted  by  this  option. 
Background  conservative  concentrations  at  each 
of  the  system  co-ordinates  are  calculated  by 
<■  linear  regression  on  simulated  water  quality 

monitoring  station  variates: 
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Cj=A^  +BjCj<  +  t*ST!).KRR.EST.i  i-coordinat©  index 

k-water  quality  monitoring 
station  index 

where  A^B- ,STD.ERR.KST. ^  are  input 
t-N(O.l) 

If  any  of  tlie  test  points  or  junction  points  are  to 
consider  excess  flows  above  a  pre-specified  target 
flov;  (channel  capacity),  then  punch  1.  The  same 
statistical  summary  output  for  excesses  or  deficiencies 
will  be  provided. 

If  evaporation  at  dam  sites  is  desired,  punch  1. 

Printing  suppress  option.  Punch  1  if  not  desired. 

Seasonal  supplemental  flow  histograms  associated 
with  pollutant  1. 

Seasonal  supplemental  flow  histograms  associated 
with  pollutant  2. 

Seasonal  supplemental  flow  histograms  associated 
with  pollutant  3. 

Seasonal  supplemental  flow  histograms  associated 
with  pollutant  4. 

Seasonal  supplemental  flow  histograms  associated 
with  pollutant  5. 

Seasonal  supplemental  flow  histograms  associated 
with  pre-specified  target  flows. 

Maximum  seasonal  supplemental  flow  histograms. 

Seasonal  means,  standard  deviations  and  probabilities 
of  supplemental  flows  associated  with  pollutant  1. 

Seasonal  means,  standard  deviations  and  probabilities 
of  supplemental  flows  associated  with  pollutant  2. 

Seasonal  means,  standard  deviations  and  probabilities 
of c supplemental  flows  associated  with  pollutant  3. 


Seasonal  means,  standard  deviations  and  probabilities 
of  supplemental  flows  associated  with  pollutant  4. 

Seasonal  means,  standard  deviations  and  probabilities 
of  supplemental  flows  associated  with  pollutant  5. 

Seasonal  means,  standard  deviations  and  probabilities 
of  supplemental  flows  associated  with  prc-spccified 
target  flows. 

Maximum  seasonal  means,  standard  deviations  and 
probabilities  of  supplemental  flows. 

Yearly  means,  standard  deviations,  probabilities 
and  histograms  of  supplemental  flows. 

Seasonal  means  and  histograms  of  reservoir  storage. 

Seasonal  means,  standard  deviations,  and  histograms 
of  system  flows.  Grand  means  and  standard  deviations 
of  system  flows. 

Seasonal  means  and  standard  deviations  of  pollutant 
concentrations  at  test  sites.  Time-invariant  histograms. 

Time-invariant  histograms  of  pollutant  concentrations 
at  dam  sites.  Seasonal  means  and  standard  deviations 
of  pollutants. 

Description  of  first  coordinate. 

Description  of  second  coordinate. 


Description  of  last  coordinate. 
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C.  Design  Cards 

Card  Columns 

*11  Up  to  three  cards  may  be  required. 

Order  of  the  dam  sites  has  already  been  established 
by  cards  *4  and  7.  Follow  that  order. 

I- 10R.  Capacity  of  dam  l(same  units  as  stream  flow). 

II- 20R.  Capacity  of  dam  2. 


41-50R.  Capacity  of  dam  5. 

1-10R.  Capacity  of  dam  6. 


4  *  *  •  1  Capacity  of  last  dam. 

*12  Up  to  three  cards  may  be  required, 

I- 10R.  Dead  storage,  dam  ]. 

II- 20R.  Dead  storage,  dam  2. 


41-50R.  Dead  storage,  dam  f 

•  •  •  • 

1-10R.  Dead  storage,  dam 


•  *  •  •  ‘  Dead  storage,  last  dam. 

*13  1-10R.  Initial  storage,  dam  1. 

11-20R.  Initial  storage,  dam  2. 


I 

I 

I 

I 

I 

l 

I 

I 

I 

I 


41-50R.  Initial  storage,  dam  5. 

•  •  •  «  • 

1-10K.  Initial  storage,  dam  6. 


•*  *e*  Initial  storage,  last  dam. 

NOTE:  Units  of  streamflows,  capacities,  diversions, 
flow  requirements,  and  waste  volumes  must  be 
made  consistent  external  to  the  program. 

Reservoir  rule  curves. 

General  expression  for  draft  during  any  season: 

=  Aj+BiQ  +  CjSj  i  =  seasonal  index 

where  A*.,  Bj,  are  input 

If  B^  $  =  0,  and  A^  less  than  0  ,  (negative 

decimal)  then  a  percentage  of  capacity(pool  level 
rule)  is  implied. 

I- 10R.  Coefficient  A,  season  1,  dam  1. 

II- 20R.  Coefficient  A,  season  2,  dam  1. 


51-60R.  Coefficient  A,  season  C.  dam  1. 


1-10R*  Coefficient .A,  season  7,  dam  1. 


51-60R.  Coefficient  A,  season  12,  dam  1. 


1-10R.  Coefficient  B,  season  1,  dam  1. 


51-60R.  ^Coefficient  B,  season  6,  dam  1. 

1 


1 - 1  OR . 


Coefficient  B,  season  7,  dam  1. 


51-6011. 


Coefficient  B,  season  12,  dam  1. 


1-10!’. 

11-20R. 

51-60R. 

1-10R. 

51-60R. 

4-5R 

9-1  OR 

49-50R 

4-5R 


\ 
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Coefficient  C,  season  1,  dam  1. 
Coefficient  C,  season  2,  dam  1. 


Coefficient  C,  season  6,  dam  1. 


Coefficient  C,  season  7,  dam  1. 


Coefficient  C,  season  12,  dam  1. 


Cycle  for  all  dans. 


If  explicit  (tvansbasin)  diversion  choice  is 
exercised  (8-55),  include  this  card(s). 

Location  of  test  site  (t~st  points  and  confluence 
points)  that  dar.t  1  is  to  make  diversions  to. 

Location  of  test  site  that  dam  2  is  to  make 
diversions  to. 


Location  of  test  site  that  dam  10  is  to  make 
diversions  to. 

Location  of  test  site  that  dam  11  is  to  make 
diversions  to. 

Location  of  test  site  that  last  dam  is  to  make 
diversion  to. 

^Explicit  mean  divers’ on  for  each  dam. 


•  •  •  « 


I- 10U. 

II- 20U. 


51-60R. 


1-10R. 


S1-60R. 


4-5R 


49-SOR 


4-5R 
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Diversion,  dam  1,  season  1. 
Diversion,  dam  1,  season  2. 


Diversion,  dam  1,  season  6. 


Diversion,  dam  1,  season  7. 


Diversion,  dam  1,  season  12. 

Cycle  for  all  dams. 

Units  must  be  made  consistent  external  to  program. 

If  conditional  diversion  option  (8-40)  is  exercised, 
include  this  card(s). 

Location  of  test  site  that  dam  1  will  supplement 
with  diversions  when  needed. 


Location  of  test  site  that  dam  10  will  supplement 
with  diversions  when  needed. 


Location  of  test  site  that  dam  11  will  supplement 
with  diversions  when  needed. 


Location  of  test  site  that  last  dam  will  supplement 
with  diversions  when  needed. 

NOTF.:  If  no  waste  variable  is  specified  (although 
waste  sites  may  be  specified)  skip  to  card  #21. 

All  waste  volumes  are  waste  flow  rates.  All  waste 
weights  arc  the  waste  volume  rates  times  the 
concentration  of  the  respective  pollutants.  Order 
of  waste  sites  has  already  been  established  by  cards 
#4  and  7.  Folic,  that  order. 


40 


18 


I- 1  OR. 

II- 20R. 


41-50R. 


Minimum  volume,  site  1. 
Minimum  volume,  site  2. 


Minimum  volume,  site  5. 


19 


*20 


I- 10R. 

II- 20R. 


41-50R. 


I- 10R. 

II- 20R- 


51-60R. 


I- 10R. 

II- 20R. 


Minimum  volume,  last  site. 

Minimum  pollutant  weight,  waste  1,  site  1. 

Minimum  pollutant  weight,  waste  2,  site  1. 


Minimum  pollutant  weight,  waste  5,  site  1. 


Start  new  waste  site. 

Note:  Cards  *18  and  *i9  constrain 
stochastic  waste  pollutant  volume  and 
weights  to  minimum  input  values.  Punch 
zeros  if  the  respective  standard  deviations 
in  *20  are  zero. 

Seasonal  moments  of  waste  volume. 

Total  mean  volume,  waste  site  1,  season  1. 
Total  mean  volume,  waste  site  1,  season  1. 


Total  mean  volume,  waste  site  1,  season  6. 

Total  mean  volume,  waste  site  1,  season  7. 

Total  mean  volume,  waste  site  1,  season  8. 


Total  mean  volume,  waste  site  1,  last  season. 


4] 

I- 10R  •  Standard  deviation  of  volume,  waste  site  1,  season  J. 

II- 20R.  Standard  deviation  of  volume,  waste  site  1,  season  2. 


51-60R.  Standard  deviation  of  volume,  waste  site  1,  season  6. 


I- 10R,  Standard  deviation  of  volume,  waste  site  1,  season  7. 

II- 20R,  Standard  deviation  of  volume,  waste  site  1,  season  8. 


Standard  deviation  of  volume,  waste  site  1,  last, 
season. 


Seasonal  moments  of  pollutant  weights. 

I- 10R.  Mean  weight  of  waste  1,  season  1. 

II- 20R.  Mean  weight  of  waste  1,  ,  season  2. 


51-60R.  Mean  weight  of  waste  1,  ,  season  6. 

1-10R.  Mean  weight  of  waste  1,  ,  season  7. 


Mean  weight  of  waste  1,  last  season. 


1~10R.  Standard  deviation  of  weight  for  waste  1,  season  1. 

11-20R.  Standard  deviation  of  weight  for  waste  1,  season  2. 


51-CdR.  Standard  deviation  of  weight  for  waste  1,  season  6. 

1-10R.  Standard  deviation  of  weight  for  waste  1,  season  7, 


fl 
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T 


.  Standard  deviation  of  weight  for  waste  1,  last  season 

Repeat  for  each  waste  variable  then  start  at  next 
waste  site. 

Three  cards  may  he  necessary. 

*21  Number  of  1st  tier  dams  upstream  servi;^ 

downs tr  ma  test  points  and  junction  points  U  10). 
Order  c f  test  sitos  has  already  been 
established  by  cards  *4  and  7.  Follow  that 
order.  * 


1-5R 

Number  of  dams  serving  test  point  or  junction 
point  1. 

6-1CR 

Number  of  dams  serving  test  point  or  function 
point  2. 

46-50R 

Number  of  dams  serving  test  point  or  junction 
point  10. 

1-5R 

Number  of  dans  serving  test  point  er  junction 
point  11. 

Number  of  dans  serving  last  test  point 
(junction  point). 

25  cards  may  be  required.  Location  cf  dams 
serving  test  points.  If  no  dans  serve  test 
point,  then  skip  that  card. 

*22  1-5R 

location  of  dam  1  serving  test  point  1. 

6-1  OR 

Location  of  dan  2  serving  test  point  l. 

11-1SR 

Location  of  dam  5  serving  ust  point  1. 

Location  of  last  dam  serving  test  point  1. 

1  -SR 

Location  of  dan  1  serving  test  point  2. 

•  »  *  • 

Location  of  last  dan  serving  last  test  point. 

•Junction  points  arc 

included  in  the  ordering  after  both  branch 

b^e!i  co  pi  etc J 

« 

i 


43 


23  35U  cards  may  bo  necessary.  Seasonal  moments 

of  flow  requirement  and  seasonal  targets  for 
each  waste  variable  at  each  test  point. 

Test  Point  1. ~ 


1-10FL 

- Mean  flow  requirement,  season  1. 

11-20R. 

Mean  flow' requirement,  season  2. 

51-60R. 

Mean  flow  requirement,  season  6. 

1-10R. 

Mean  flow  requirement,  season  7. 

•  >  ■»  t  * 

Mean  flow  requirement,  last  season. 

l'-ibft. 

Standard  deviation  flow  requirement , 

season  1. 

11-20R. 

Standard  deviation  flow  requirement. 

season  2. 

S1-60R. 

Standard  deviation  flow  requirement. 

season  6. 

1-10R. 

Standard  deviation  flow  requirement. 

season  7. 

•  « 

Standard  deviation  flow  requirement. 

la: t  ^eas^n. 

Water  quality  requirements  (Standards)  for  carh 
test  point.  Units  must  be  consistent  with  pollutant 
concentrations.  Note;  Skip  target-cards  if  no 
pollutant  specified. 

1-10R. 

Target  waste  1,  season  1. 

11-20R. 

Target  waste  !,  season  2. 

51-60R. 

Target  waste  1,  season  6. 

1-10R. 

Target  waste  1,  season  7, 

Target  waste  1,  last  season. 

Finish  all  wastes,  then  start  the  sane 
sequence  st  next  test  point  on  junction 
point. 

Option  "to  keep  track  of  flows  higher  than 
a  specified  liait.  (Flood  flows) 

If  8-60  was  non-iero,  then  include  *24 

If  the  flow  requirement  at  test  point  or  junction  - 
point  1  is  a  channel  capacity  restraint,  punch  1. 

If  the  flow  requirement  at  test  point  or  junction 
point  2  is  a  channel  capacity  restraint,  punch  1. 


If  the  flow  requirement  at  test  point  or  junction 
point  10  is  a  channel  capacity  restraint,  punch  1. 

Evaporation  pararctar  input. 

If  card  8-65  was  zero  then  skip  to  *53. 

2  cards  cay  be  necessary. 

Punch  1  if  contents  of  dam  1  are  to  be  evaporated. 
Punch  1  if  contents  of  das  2  are  to  be  evaporated. 


Punch  1  if  contents  of  dan;  10  are  to  be  evaporated.. 

Punch  1  if  contents  of  last  dam  arc  to  be  evaporated 

Number  of  regional  seasonal  evaporation  vectors,  £3, 
(No  conversion  of  units  is  performed  in  progran). 


I 

I 

I 

e 

o 


I- 10R.  Mean  evaporation  rate  for  season  1,  vector  1. 

II- 20R.  Mean  evapor.  tion  rate  for  season  2,  vector  1. 


51-60R.  Mean  evaporation  rate  for  season  6,  vector  1. 


.  Mean  evaporation  rate  for  last  season,  vector  1. 

I- 1CR.  Standard  deviation  of  evaporation  rate  for  season  1 , 

vector  1. 

II- _0R.  Standard  deviation  of  evaporation  rate  for  season  2, 

vector  1. 


S1-60R.  Standard  deviation  of  evaporation  rate  for 

season  6,  vector  1. 


Standard  deviation  of  evaporation  rate  for  iast 
season,  vector  1. 

Repeat  cards  *27  and  *28  for  next  evaporation  vector. 


*29  Two  cards  may  be  needed. 

5  Evaporation  vectors  (mean  and  standard  devi<tlion) 

to  which  the  first  dam  is  to  be  referred. 

10  Evaporation  vector  to  which  the  second  dam  is  to 

be  referred.  ’ 


SO 


Evaporation  vector  to  which  10th  dam  is  to  be 
referred. 


Evaporation  vector  to  which  last  dam  is  to  be 
referred. 
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I- 10R.  Initial  concentration  of  waste  1,  1st  site 

with  storage?©, 

II- 20R.  Initial  concentration  of  waste  2,  1st  site 

with  storage y  0. 

21-30R.  Initial  concentration  of  waste  3,  1st  site 

with  storage  7 0. 

31-40R.  Initial  concentration  of  waste  4,  1st  site 

with  storage  7  0. 

41-50R.  Initial  concentration  of  waste  5,  1st  site 

with  storage y  0. 

1-10R.  Initial  concentration  of-vaste  1,  2nd  site 

with  storage?©. 


41-SOR.  Initial  concentration  of  waste  5,  last  site 

with  storage  ?0. 

5  If  waste  1  is  degradable  punch  1;  if  not,  zero 

10  Number  of  conservative  wastes. 

a 

If  1  is  punched  in  column  5  of  card  34,  the 
following  cards  are  required.  If  zero,  skip 
to  *42. 


As  many  as  9  cards  may  ae  required. 

I- 10R.  Non-conservative  pollutant  concentration  of 

unregulated  flow,  point  1. 

II- 20R.  Non-conservative  pollutant  concentration  of 

unregulated  flow,  point  2. 


51-60R.  Non-conservative  pc’Hitant  concentration  of 

unregulated  flow,  point  6. 

•  •  -  «  • 

1-10R.  Non-conscrvativc  pollutant  concentration  of 

unregulated  flow,  point  7. 


I- 1 OR. 

II- 20R. 


51-60R. 

i-iop.. 


I- 1 OR. 

II- 20R. 


51-60R. 
•  *  \  ' 

1-lOR. 


•  •  •  • 


I- lOR. 

II- 20R, 


51-60R. 

1-lOR. 


Non-conservative  pollutant  concentration  of 
unregulated  flow,  last  point. 

Two  cards  may  be  required. 

Mean  temperature  (°C),  season  1. 

Mean  temperature  (°C),  season  2. 


Mean  temperature  (°C) ,  season, 6. 

Mean  temperature  (°C),  season  7. 

Mean  temperature  (°C),  last  season. 

Nine  cards  may  be  required. 

River  mile  ..e  from  point  to  next  point. 
River  mileage  from  point  2  to  next  point. 


River  mileage  from  point  6  to  next  point. 
River  mileage  from  point  7  to  next  point. 

River  mileage  last  point,  punch  zero. 


Nine  cards  may  be  required.  Stream  velocity 
is  computed  from:  V=aQ^. 

Constant  "a"  in  equation  for  velocity  at  site  1. 
Constant  "a"  in  equation  for  velocity  at  site  2. 


Constant  "a"  in  equation  for  velocity  at  site  6. 
Constant  "a"  in  equation  for  velocity  at  site  7. 


in  equation  for  velocity  at  last  site 


Constant  "a" 


I- 10R. 

II- 20R. 


51-60R. 

1-10R. 

I- 10R. 

II- 20R. 


I- 10R. 

II- 20R. 


51-60R. 
1  - 1  OR . 


Nine  cards  nay  be  required. 

Constant  "b"  in  equation  for  velocity  at  point  1 . 
Constant  "b"  in  equation  for  velocity  at  point  2. 


Constant  "b"  in  couation  for  velocity  at  point  6. 
Constant  "b"  in  equation  for  velocity  at  point  7. 


Deoxygenation  coefficient  (20°C)  at  point  1. 
Ceoxyger.ation  coefficient  (206C)  at  point  2. 


Deoxygenation  coefficient  (20cC)  at  last  point. 
Three  cards  nay  be  required. 

Deoxygenation  coefficient  (20°C)  at  dam  site  1. 
Deoxygenation  coefficient  (20°C)  at  dam  site  2. 


Deoxygenation  coefficient  (20°C)  at  dam  site  6. 
Deoxygenation  coefficient  (20°C)  at  dam  site  7. 


If  34-10  is  zero,  skip  to  *45 
Conservative  pollutant  background  concentrations. 
If  S-55  was  punched  1,  use  Option  1:  if  2,  use 
Option  1  option  2-  The  expected  value  of  the  concen¬ 
tration  for  each  conservative  pollutant  at  all 
sites  is  calculated  from  the  exponential 
regression  on  flow: 

!_  '  C=dQf  C  T 


Ole 


t  *  mawmnr  iwwMKMjr 


...  X 
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Option  1 

I- 10R.  "d"  in  the  equation  for  concentration  1,  site  1. 

II- 20R.  "d"  in  the  equation  for  concentration  2,  site  1. 

•  «•••  ••••• 

41-50K.  "d"  in  the  equation  for  concentration  5,  site  1. 

I- 1  OR.  "f"  in  the  equation  for  concentration  1,  site  1. 

II- 20R.  "f"  in  the  equation  for  concentration  2,  site  1. 

41-50R.  "f"  in  the  equation  for  concentration  5,  site  1. 


I- 10P.  Maximum  value  ("T")  of  concentration  1,  site  1. 

II- 20P..  Maximum  value  ("T”)  of  concentration  2,  site  1. 


41-50R.  Maximum  value  ("T")  of  concentration  5,  site  1. 

I- 10R.  Minimum  value  ("L”)  of  concentration  1,  site  1. 

II- 20R.  Minimum  value  ("L")  of  concentration  2,  site  1. 


41-50R.  Minimum  value  ("L")  of  concentration  5,  site  1. 

Cycle  for  all  sites 

Option  2  Option  limited  to  one  conservative  pollutant. 

Simulated  stream  flows  and  water  quality  are  on 
the  flow  (quality)  tape.  Caution:  check  #5  for 
consistency  with  *43.  All  entries  in  *43  must 
be  different  than  those  in  *5. 

*43  4-5R  Water  quality  monitoring  station  which  coordinate 

1  is  referred. 

9-10R  Water  quality  monitoring  station  which  coordinate 

2  is  referred. 
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14-13R 


49-50R 


4-5R 


*44 

1-10R. 


1-10R. 


1-10R, 


*45  1-6R. 


7-12R, 


\ 
i 

\ 

i 

i 

i 

i 

Water  quality  monitoring  station  which  coordinate! 

3  is  referred. 


Water  quality  monitoring  station  wich  coordinate 
10  is  referred. 


Water  quality  monitoring  station  which  coordinate 
11  is  referred. 

Water  quality  monitoring  station  which  last 
coordinate  is  referred. 

Simulated  values  of  the  conservative  pollutant  are 
calculated  for  all  sites  by  linear  regression  on 
water  quality  monitoring  station  variates.  See-8-55 

Constant  in  equation  for  concentration  1  at.  site  1. 


Regression  coefficient  in  equation  for  concentration  j 
1  at  site  1.  I 


Standard  Error  of  Estimate  in  equation  for  concentra¬ 
tion  1  at  site  1. 

Cycle  for  all  sites. 

Up  to  five  cards  may  bo  required. 

Initial  flow  at  gaging  station  1. 

Initial  flow  at  gaging  station  2. 


i 
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If  input  hydrology  is  on  cards, 

flows  at  gaging 

stations  in  chronological  order 

Year  1  Season  1. 

are  needed. 

1-6R. 

Flow  at  gaging  station  1. 

7-12R. 

Flow  at  gaging  station  2. 

67-72R. 

Flow  at  gaging  station  12. 

1-6R. 

Flow  at  gaging  station  13. 

••••»  Flow  at  last  gaging  station. 


last  year  final  season. 
1-6R.  Flow  at  gaging  station  1. 

7-12R.  Flow  at  gaging  station  2. 


67-72R.  Flow  at  gaging  station  12. 


.  Flow  at  last  gaging  station. 

*47  Data  Supplement  Option  (Multi-run  specification) 

Comment  cards.  Five  are  needed, 

1-72  Comment  card  1. 


Comment  card  5. 
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*48 


48-A  4~5R 

9-10R 


14-15R 


19-20R 


24-25R 

29-30R 

34-35R 

39-40R 

44-45R 

49-50R 


Two  cards  are  necessary.  Number  of  substitutions 
for  each  of  the  "design"  and  "data"  parameters. 

No  geometry  card  substitutions  arc  permitted. 

A  complete  seasonal  replacement  is  counted  as 
one  substitution. 

Design  cards. 

Number  of  reservoir  volume  changes. 

Number  of  dead  storage  changes. 

Number  of  starting  storage  changes. 

Reservoir  rule  curve  change. 

Draft  -  A  ♦  B*Q  +  C*S 

Number  of  changes  of  the  seasonal  coefficient  "A" 

Note:  Complete  seasonal  change  is  counted  as  one 
substitution.  Pool  level  rules  are 
negative  fractions. 

Number  of  changes  of  coefficient  "B"  for  all 
reservoirs . 

Number  of  changes  of  coefficients  "CV  for  all 
reservoirs. 

Number  of  test  points  and  junction  points  for  which 
candidate  dams  for  augmentation  are  to  be  changed. 

Number  of  seasonal  explicit  (transbasin)  diversions 
to  be  changed. 

Number  of  test  and  junction  points  that  ore  to  receive 
conditional  diversions  for  a  different  dam. 

Number  of  waste  volume  means  to  be  changed.  A  complete 
seasonal  change  at  the j  1  waste  site  is  counted  as 
one  change. 


48-B  4-5U  Number  of  waste  volume  standard  deviations  to  bo 

changed. 


9-1  OR 


Number  of  mean  pollutant  weights  to  be  changed. 


i 


*49-A 


S4 


14-15R  Number  of  pollutant  weight  standard  deviations 

to  be  changed.  (A  complete  seasonal  change  for  the 
ith  pollutant  at  the  jtil  waste  site  is  counted  as 
one  change) . 

19-20R  Number  of  mean  low  flow  requirements  to  be  changed. 

(A  complete  seasonal  change  at  the  j*'1  test  site 
is  counted  as  one  change) . 

24-25R  Number  of  flow  requirement  standard  deviations 

to  be  changed. 

29-30R  Number  of  water  quality  target  changes.  (A  complete 

seasonal  vector  change  for  the  it  1  pollutant  at  the 
jth  test  site  is  counted  as  one  substitution. 

Data  Cards 

Conservative  pollutant  background. 

Option  A:  L ^  C  -  T 

Option  E:  C  «  d  ♦  fC  ♦  S.E.E.  *t  t  -  (0,1  )N 
34-38R  Number  of  constant  "d"  changes.  (Option  l  or  2). 

39-40R  Number  of  exponent  "f"  changes.  (Option  1  or  2). 

44-4SR  Number  of  maximum  value  "T"  changes.  (Option  1  only) 

Number  of  "S.F..F.."  changes.  (Option  2  only). 

49-50R  Number  of  minimum  value  ”L"  changes.  (Option  l  only) 

Dam  capacity  Changes.  (Refer. to 
4-SR  Daw  number  (not  coordinate). 

6-15U.  Dam  capacity. 

Repeat  for  as  many  changes  indicated  by  48-A-5 
Dtad  r  tor  chany.es .  (Refer  to  *48- A) 

4-5R  Daa  number  (not  coordinate). 

6-15R.  Deed  storage. 

Repeal  for  as  many  changes  indicated  by  M8-A-10. 
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4- 5R 

5- 15R. 


4-5H 

9-10R 

11-20R. 


I- 10R. 

II - 20R. 


51-60R. 


1-10R. 


4-5R 
9-1  OR 


Initial  storage  changes.  (Refer  to  M8-A) 

Dam  number  (not  coordinate) . 

Initial  storage. 

Repeat  for  as  many  changes  indicated  by  48-A-1S. 

Reservoir  rule  curve  coefficient  changes.  Refer  to 
to  *4 8 -A) 

Dam  number 

Season  number  (If  entire  vector  is  to  be 
replaced,  punch  number  of  seasons  ♦  1.  Then 
ship  11-20R. ) . 

Constant  "A  .  Pool  level  rules  are  negative 
fractions. 

If  change  desired  for  all  seasons,  two  cards 
may  be  necessary. 

Constant  "A",  season  1. 

Constant  A",  season  2. 


Constant  '’A",  Season  6. 

Constant  "A",  Season  7, 


Constant  "A",  last  season. 

Cycle  for  as  many  changes  indicated  by  M8-A-20. 
Dam  numbor 

Season  number  (If  entire  vector  is  to  be  rcolnced 
punch  numbor  of  seasons  ♦  1.  Then  skip  11*20!*, 


11-20R. 


*  *  '.  t  < 


1-10R. 


51-60R. 


1-10R. 


4-5R 

9-10R 

11-20R. 

1-10R. 


►1  P«J 


Regression  coefficient  "B". 

If  change  desired  for  all  seasons,  two  cards 
may  be  necessary. 

Coefficient  "B",  season  1. 


Coefficient  "B",  season  6, 

Coefficient  "B"  season  7. 

Coefficient  "B",  last  season. 

Cycle  for  as  many  changes  indicated  by  *48-A~25. 
Dam  number. 

Season  number. 

Regression  coefficient  "C",  skip  if  changing 
entire  vector. 

If  replacing  entire  seasonal  vector,  two 
cards  may  be  necessary. 

Regression  coefficient  "C",  season  1. 


51-60R, 


Regression  coefficient  "C",  season  6. 


1-10R. 


Regression  coefficient  "C", season  7. 

Regression  coefficient  "C",  last  season. 

Cycle  for  as  many  changes  indicated  by  *4$-A-30. 
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4-5R 

9-1QR 

4-5R 

9-10R 


4-5R 

9-1QR 

11-15R 

16-25R. 

•  *  5  •  4 

I- 10R. 

II- 20R. 

S1-60R. 

1-10R. 


Pre-specified  reservoirs-  (refer  to  48-A) 

Test  point  or  junction  point  number  (not 
coordinate) . 

Numoer  of  first  tier  ( 6 10)  dams  serving 
this  point. 

Location  of  1st  dam  serving  this  point. 
Location  of  2nd  dam  serving  this  point. 


Repeat  for  as  many  changes  indicated  by 
•48-A-35 

Explicit  diversion  changes.  (Refer  to  *48-A) 

Dam  number 

Location  of  test  site  that  dam  is  to  make 
diversions  to. 

Season  number  (If  replacing  entire  vector,  skip 
16-25R.) . 

Seasonal  diversion. 

If  replacing  the  entire  vector,  two  cards 
may  be  necessary. 

Diversion,  season  1. 

Diversion,  season  2. 

Diversion,  season  6. 

Diversion,  season  7- 


Diversion,  last  season. 

Cycle  for  as  many  changes  indicated  by 
M5-A-40. 
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4-5R 

9-10R 


4-SR 

9-10R 

1J.-20R. 


1-10R. 


51-60R. 


1-10R. 


Conditional  diversion  changes,  (refer  to  *48-A) . 
Dam  number. 

•w 

Location  of  test  site  receiving  supplementation 
as  needed. 

Repeat  for  as  many  times  as  indicated  in 
*48-A-45. 

Pollutant  volume  and  weight  changes! 

Mean  waste  volume  changes  .(refer  to  *48-A) 

Waste  site  (not  coordinate). 

Season  number.  (If  entire  seasonal  vector 
is  to  be  replaced,  punch  number  of  season  +1. 
Then  skip  11-20R.) 

Mean  waste  volume  change. 


If  change  desired  for  all  seasons,  then  two 
cards  may  be  required. 

Mean  waste  volume,  season  1. 


Mean  waste  volume,  season  6. 


Mean  waste  volume,  season  7. 


Mean  waste  volume,  last  season. 


Rcpcnst  for  as  many  changes  indicated  by 
*48-A-50. 


M9-B 


4-5R 


Waste  volume  standard  deviation  changes. 
Waste  site  (not  coordinate). 
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Q 

0 

D 


9-10R 

11-20R. 


Season  number.-  (if  entire  vector  is  to  bo 
replaced,  punch  number  of  seasons  +1.  Then 
skip  11-2011). 

Standard  deviation  waste  volume  change. 


If  change  desired  for  all  seasons,  then 
two  cards  may  be  necessary. 

1-10R.  Standard  deviation  waste  volume,  season  1. 


•  •  •  •  •  •  •  »  i 

51-60h.  Standard  deviation  waste  volume,  season  6. 


1-10R, 


Standard  deviation  waste  volume,  season  7. 


4-5R 
9-1  OR 


15 

16-25R. 


Standard  deviation  waste  volume,  last  season. 

Cycle  for  as  many  changes  indicated  by  *48-5-5. 

Mean  waste  weight  changes.  (Refer  to  *4S-B) 

Haste  site  (not  coordinate). 

Season  number,  (If  the  entire  vector  is  to  be 
replaced,  punch  number  of  seasons  +1,  Then 
Skip  if  16-25R, . 

Pollutant  number. 

Mean  waste  weight  change. 


If  change  desired  for  nil  seasons,  then  two 
cards  may  bo  required. 

1-10R.  Mean  waste  weight,  season  1. 


5J-60R. 


Mean  waste  vw,.  ,  season  6. 


'  31 


n 


1-1 OK . 


4-5R 
9-1  OR 

16-25R, 


I- lOR. 

II- 20R. 


51-60R. 


1-lOR. 


4-5R 
9-1  OR 

Ii-20R. 


Mean  waste  weight,  season  7. 


Mean  waste  weight,  last  season. 

Cycle  for  as  many  changes  indicated  by 
M8-B-10. 

Waste  weight  standard  deviation  changes. 

Waste  site  (not  coordinate). 

Season  number.  (If  entire  vector  is  to  be 
replaced,  punch  number  of  seasons  +1.  Then 
skip  J6-25R.) . 

Standard  deviation  waste  weight. 

If  change  desired  for  all  seasons,  two  cards 
may  be  necessary. 

Standard  deviation  waste  weight,  season  1. 
Standard  deviation  waste  weight,  season  2. 


Standard  deviation  waste  weight,  season  6. 


Standard  deviation  waste  weight,  season  7. 


Standard  deviation  waste  woltht,  last  season. 

Cycle  for  as  many  changes  indicated  by  *4S-B-15. 

Mean  flow  requirement  changes. (Refer  to  MS-R) 

Test  point  or  junction  point  (not  coordinate). 

Season  number.  (If  all  seasons  are  to  be  changed 
then  punch  number  of  seasons  *1).  Skip  next 
field.) 

Mean  target  flow. 


H 


61 


1-lOR. 

51-60R. 

1-lOR. 


4-5R. 

9-10R 

ll-JOP. 

1-lOR. 

51-60R. 

1-lOR. 


If  change  Jesired  for  all  reasons,  two 
cards  may  bo  necessary. 

Mean  target  flow,  season  1. 


Mean  target  flow,  season  6. 


Mean  target  flow,  season  7, 


Mean  target  flow,  last  season. 

Cycle  for  as  many  changes  indicated  by 
♦48-B-20. 


Flow  requirement  standard  deviation  chances. 

(Refer ‘to  *48-65 - 

Test  point  or  junction  point  (not  coordinate). 

Season  number.  (If  all  seasons  are  to  be 
changed,  then  punch  number  of  seasons  +1. 
Skip  next  field.) 

Standard  deviation  of  flow  target. 


If  change  desired  for  all  soasons,  two  cards 
may  be  necessary. 

Standard  deviation  flow  target,  soason  I. 


Standard  deviation  flow  target,  soason  6. 


Standard  deviation  flow  target,  season  7. 


Standard  deviation  from  flow  target,  last  season. 
CycU  for  as  many  changes  indicated  by  *4«-fc-2S. 
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4-5K 

9-1011 

14-15R 

16-25R. 


1-10R. 

U-20R. 

SI-60R. 

i-ioR. 


Option  ^ 


4-SR 

9-10R 

1J-20R. 


Water  quality  requirement  changes. (Refer  to  M8-B) 
Test  point  or  junction  point  (not  coordinate). 
Pollutant  number. 

Season  number.  (If  entire  seasonal  vector 
is  to  be  replaced,  punch  number  cf  seasons  +1. 
Skip  next  field). 

Value  of  target  concentration  to  be  changed. 

If  change  desired  for  all  seasons,  then  two 
cards  may  be  necessary. 

Target  concentration,  season  1. 

Target  concentration,  season  2. 

Target  concentration,  season  6. 

Target  concentration,  season 


Target  concentration,  last  season. 

Cycle  for  as  many  changes  indicated  by  M8-B-30. 
Conservative  pollutant  background  changes 

Constant  in  the  exponential  regression 
equation.  (Refer  to  *48*61 

Coordinate 

Pollutant  number. 

Constant  Md*'  change. 

Cycle  for  as  many  charges  as  indicated  by 
♦48-B-35. 


Regression  coefficient. 

Coordinate 
Pollutant  number. 

Exponent  "fM. 

Cycle  for  as  many  changes  as  indicated  by 
•48-B40. 

High  concentration  constraint  changes. 
Coordinate. 

Pollutant  number. 

Constraint  VT". 

Cycle  for  as  many  changes  indicated  by 
•48-B-4* 

Lower  concentration  constraint  changes. 
Coordinate. 

Pollutant  number. 

Constraint  “l". 

Cycle  for  as  many  changes  indicated  by  M8-B-S0. 


Option  2 


4-5R 

9-10R 

11-20R. 


4-5R 

9-10R 

11-2QR. 


4-5R 

9-lOR 

11-20R. 


Constant  in  the  linear  regression  equation 
Coord inate 

Pollutant  number  (limited  to  one) 

Constant  ”d''  change. 

Cycle  for  as  many  changes  indicated  by  *4S-B~35. 

Regression  Coefficient 
Coordinate 

Pollutant  number  (limited  to  one) 

Regression  coefficient  "f” 

Cycle  for  as  many  changes  as  indicated  by 
•48-B-40. 

Standard  Error  of  Estimate 
Coordinate 

Pollutant  number  (limited  to  <re) 

Standard  Error  of  Estimate 

Cycle  for  as  many  changes  as  indicated  by 
*48-9-45. 


GENERAL  NOTES  CONCERNING  USAGE  OF  OPTION  D 


1.  The  dam  site  and  the  use  point  should  be  proximate  to  the  gaging 
station;  however,  if  the  basin  is  very  homogeneous  with  respect  to  rainfall, 
evaporation,  and  morphology,  then  this  restriction  can  be  relaxed. 

2.  The  gaging  station  need  not  be  downstream  or  even  on  the  same 
nainstem  or  tributary  as  the  dam  site  and/or  the  use  point.  All  three 
may  be  located  on  three  different  tributaries  or  mainstems. 

3.  The  flow  at  the  gaging  station  should  be  arrayed  so  that  the  low 
flow  seasons  are  in  the  middle  of  the  year  while  the  high  flow  seasons 
begin  and  end  the  year.  This  facilitates  computing  annual  supplemental 
storage  requirements. 

4.  All  units  of  flows,  storage  and  requirements  must  be  made 
compatible,  external  to  the  program. 

5.  Perform  a  routing  using  a  reservoir  with  capacity  equal  to 
zero  before  selecting  a  number  of  trial  reservoir  sizes  to  meet  some 
objective.  The  annual  supplemental  requirements  will  serve  as  a  guide 
in  reducing  the  number  of  trial  reservoirs.  Furthermore,  if  the  success 
probability  is  very  high  using  zero  reservoir  sire,  say  96?;,  then  the 
necessary  amount  of  storage  to  insure  higher  degrees  of  protection 

can  be  computed  directly  from  the  distributions  of  annual  supplemental 
flows.  Thus  there  would  be  no  need  for  further  routing  analysis. 

6.  For  eastern  streams  where  the  required  storage  may  be  between 
C,1  -  0.5  of  the  annual  flow  at  the  reservoir  site  and  where  the  record 
is’ abundant,  the  storage,  in  the  reservoir  should  be  stai'tcd  at  901  full. 

For  studies  in  which  carry  over  storage  is  anticipated,  the  sensitivity 
of  the  starting  storage  should  be  parametrically  considered. 

7.  If  the  average  of  the  flow  requirements  is  equal  to  or  greater 
than  the  mean  flow  at  the  test  site,  then  expect  at  least  501  failures 
using  any  trial  reservoir  size. 

fc.  If  a  reservoir  is  to  provide  a  high  degree  of  protection  or  a 
high  yet urn  period  (9$v  -  1/50  failure)  then  caution  should  be  exercised 
in  selecting  the  proper  length  of  tec or df or  routing.  The  length  of 
record  should  be  grcater  *han  the  drsired  return, period. 

9.  To  test  the' hydrologic  sensitivity  of; a  design  where  the  objective 
is  to  find  the  amount  of  storage  to  meet  flow  requ i resent for  a  given 
protection  level,  employ  *  no  following  strategy: 


a)  Using  the  historical  record  and  trial  reservoir  sizes,  perform 
routings  to  find  the  required  storage.  Call  this  the  central 
design. 

b)  Select  a  size  different  from  the  central  design  value.  Perform 
at  least  twenty  routings  using  different  operational  (synthetic) 
traces  for  the  same  reservoir  size  and  dilution  requirements. 

c)  Tabulate  and  plot  percentage  of  failure  vs  storage  for  each 
routing. 

d)  Select  another  size  and  perform  at  least  twenty  routings  using 
the  same  set  of  operational  traces.  Repeat  c. 

c)  At  the  desired  protection  level,  pick  off  the  twenty  storages 
and  compute  mean  and  standard  deviation. 

f)  Realistic  probability  statements  can  now  be  made  concerning  safety 
factors  or  additona?  ‘corage  to  bo  added  to  the  central  design. 
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GENERAL  NOTES  CONCERNING  THE  USAGE  OP  OPTION  NMP 


The  grand  strategy  of  gaging  station  selection  is  to  capture  the 
pertinent  hydrologic  characteristics  of  the  basin  using  the  minimal 
number  of  gages.  No  clear-cut  rules  can  be  given;  however,  the  following 
guides  may  prove  helpful : 

1.  The  gages  should  define  ns  much  of  the  area  as  possible.  Studies 
to  date  indicate  that  between  201  to  60%  of  the  area  is  sufficient. 

2.  The  gages  should  have  concurrent  record  lengths  of  at  least 

20  years.  A  sufficient  number  of  wet  and  arid  periods  should  be  included. 

3.  They  should  adequately  represent  the  regional  differences;  that  is, 
the  sub-basins  exhibiting  different  types  or  degree  of  rur.off  should  be 
included, 

4.  Gages  that  are  proximate  to  important  locations  (new  reservoirs 
or  control  points)  should  be  included.  Error  is  introduced  wherever 
information  (flows  at  gaging  station)  is  translated  or  sealed  to  define 
the  inflows  at  other  points.  Considerable  error  may  6c  introduced  when 
defining  the  inflows  at  a  reservoir  located  ir.  a  headwater  reach  from  a 
downstream  gage  whose  variability  and  correlative  s true*  ire  is  different. 

5.  Gages  that  have  been  altered  by  man-made  activities  are  not 
desirable.  Cages  that  are  downstream  of  existing  reservoirs  that  have 
detention  times  of  greater  than  R.l  cf  a  year  should  be  avoided.  If  it 
is  necessary  to  include  the  gage,  records  can  be  back-routed  to  natural 
conditions:  A  hydrologist  should  be  consulted  if  this  problem  arises. 

6.  The  inclusion  of  many  mainstem  gages  should  bo  avoided  because 
the  information  contained  in  the  record  of  tie  upstream  gages  is  also 
contained  in  the  downstream  gages.  Tins  is  redundant  information.  Side 
gages  should  be  used;  each  generally  represents  new  or  frcs'i  information. 
However,  the  shape  of  the  basin  and  length  or  concurrent  records  will 
usually  limit  this  consideration. 

7.  This  option  can  also  be  used  to  simulate  together  with  stream 
flows,  traces  of  conservative  pollutants  such  as  total  dissolved  solids, 
sulfates,  and  cMorid-s.  In  genial ,  maMts*.  m  water  quality  gages 
downstream  should  be  avoided.  The  implications  of  the  other  guidelines 
are  also  applicable. 


GENERAL  NOTES  CONCERNING  THE  USAGE  OF  OPTION  B 


1.  A  nap  defining  the  system  end  all  its  components  should  bo  drawn 
and  should  be  referred  to  at  all  tines.  This  map  should  also  delineate 
the  different  areas  of  hydrologic  and  water  quality  characteristics.  This 
will  prove  most  useful  when  modeling  and  verifying  the  natural  system. 

2.  The  system  configuration  should  be  kept  as  simple  and  ns  small  as 
possible.  In  fact,  if  it  is  possible,  complex  problems  should  be  subdivided 
into  smaller  segments.  Over  the  course  of  an  effort,  a  considerable  number 
of  man-days  will  be  saved  that  would  otherwise  be  spent  purging  errors 

from  input  decks  for  complex  systems.  Given  that  the  basic  data  such  as 
drainage  area  ratios,  target  flows,  v.'aste  loadings,  reservoir  operating 
rules  and  capacities  and  so  on  have  been  compiled,  and  that  familiarity 
with  the  option  already  exists,  new  input  decks  can  'a  re-constituted  in 
a  very  short  period  cf  time. 

3.  It  is  most  desirable  to  have  a  well-developed  plan  of  attack  for 
sorting  through  various  alternative  designs.  This  is  not  always  possible 
because  the  results  .v.  a  sometimes  not  predictable;  but  experience  has  shown 
that  a  systematic  and  orderly  mode  of  attack  will  save  time  and  money. 

Building  _thc_  hydrologic  _moclel_ 

4.  The  overall  strategy  is  twofold:  First,  the  natural  regime  is 
approximated  by  defining  the  pertinent  characteristics  at  each  of  the  discre 
points  or  locations  in  the  system.  Second,  all  significant  existing  rc-guUt* 
is  t'a m superimposed  upon  it.  Every  attempt  should  be  made  to  separate 

the  effects  of  man. 


In  modeling  the  natural  system,  several  factors  should  be  considered. 

First  of  all,  the  gccmorpholo.,)  of  the  basin  should  be  carefully  studied 
and  differences  noted.  For  example,  runoff  from  an  area  dotted  with  many 
small  lakes  is  apt  to  have  less  variability  and  have  a  stronger  temporal 
correlative  structure  than  runoff  from  terrain  of  little  relief.  Secondly, 
weather  patt'rns  and  elevation  differences  arc  also  important  and  should 
be  noted.  Yields  f : era  different  portions  ef  a  ?>•  'ugly  homogeneous  area 
can  vary  significantly.  This  information  ia  useful  in  defining  the  runoff 
From  unpaged  areas  by  using  data  from  other  portions  of  the  basin. 


A  suggested  technique  for  quickly  determining  a  suitable  Metnl  of  gag's 
is  the  following:  1}  f*v  pute  grand  roans  of  upstrea'*  and  side  gages;  2) 
Compute  yield  for  u.sgagcd  aieas  using  co  Mn.nt ions  of  the  gaged  means.  Here 
the  weighting  factors  can  be  computed  bv  drain,  q.v  area,  ratios  of,  if  side 
information  is  avail.’!  le  such  as  a  few  years  of  gaged  data,  th.cn  ratios  of 


n  flows  can  Le  r : 


tc 

with 

varii.ee. 


3)  iTYCcvd  downstream  and  cn  pare  summation  of  up'tv 


u  stem  mean;  ■!)  Recycle  am!  adjust  for  seasonal  means  and 


If  existing  reservoirs  alter  the  nature  of  the  inflows  to  them,  then 
the  reservoirs  and  their  mode  of  operation  must  be  superimposed  on  the 
natural  regime.  Reservoirs  that  have  a  detention  time  of  0.1  year  or  less 
can  generally  be  eliminated.  A  convenient  way  to  determine  whether  an 
existing  reservoir  changes  the  distribution  of  inflows  is  to  plot  the 
observed  releases  as  a  function  of  the  sum  of  the  observed  inflows  and 
beginning  of  season  storage  volumes.  If  the  line  of  best  fit  departs 
significantly  from  the  45  degree  equivalence  lino,  then  the  reservoir 
must  be  included  in  the  analysis.  The  coefficients  (seasonal)  that  define 
the  line(s)  of  best  fit  can  then  be  inputed  to  the  program  as  data.  A 
horizontal  line  implies  that  the  rule  is  to  always  maintain  a  certain  flow 
requirement.  A  line  having  slope  equal  to  one  and  displaced  from  the  origin 
implies  a  pool  elevation  rule.  The  coefficient  for  this  method  of  operation 
is  computed  by  dividing  the  horizontal  distance  (quantity)  between  the 
equivalence  line  and  the  line  of  best  fit  by  the  capacity  of  the  reservoir. 

Existing  diversions  that  are  returned  to  the  system  during  different 
portions  of  the  year  may  be  approximated  by  linking  a  test  site  with  a  waste 
site.  Irrigation  return  flows  and  industrial  and  municipal  wastes  that  are 
detained  in  lagoons  may  be  treated  in  this  fashion.  The  diversions  and 
their  return  flows  may  be  at  different  locations  in  the  basin. 


5. Drainage  area  scaling  factors.  Inflows  to  coordinates  that  are  not 
at  gaging  stations  may  be  scaled  by  a  set  of  linear  factors.  These  can  be 
computed  by  the  ratio  of  drainage  areas  or  by  the  ratio  of  mean  flows. 

It  should  be  kept  in  mind  that  these  coefficients  are  fixed  for  the  entire 
year  and  have  no  random  fluctuations  added.  The  geomorphology  and  rainfall 
patterns  are  important  in  this  respect  and  should  be  carefully  considered. 
Furthermore,  as  the  length  of  the  timeframe  decreases  the  coefficients 
of  variation  for  these  factors  will  generally  increase.  Thus,  considerable 
error  might  be  introduced  into  the  scaled  flows.  Therefore,  one  must 
decide  beforehand  which  seasons  or  periods  of  tine  arc  crucial  in  the 
analysis  and  adjust  for  them  accordingly. 

6.  Verifying  the  Model.  For  verifying  the  model  the  following 
procedure”*!?"  suggested: 

1)  Downstream  gaging  stations  should  be  included  as  test  sites  in  the 
model . 

2)  A  concurrent  period  of  record  for  both  upstream  and  downstream 
gaging  stations  must  be  selected.  If  existing  reservoirs  and  their  mode 
of  operation  have  to  be  considered,  then  the  period  of  record  should  start 
when  the  reservoirs  wore  first  operated. 

3)  Using  upstream  data,  route  the  period  of  record  down  through  tne 
system  to  produce  data  at  the  test  sites. 

4)  The  seasonal  statistics  of  the  measured  or  observed  data  at  down¬ 
stream  gaging  stations  should  be  compared  with  the  statistics  of  the 
downstream  data  produced  by  the  model.  Ihe  closeness  of  fit  chnr-*ctcrlrcs 
the  validity  of  the  model.  Tuc  difference  betvvcn  the  seasonal  r.v.nns  should 
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be  quite  small.  Often  the  seasonal  standard  deviations  produced  by 
the  model  are  less  than  the  observed  values.  This  usually  occurs 
when  a  substantial  amount  of  the  drainage  area  is  not  gaged  and  the 
runoff  has  been  defined  or- computed  from  other  gages. 

7.  Length  of  timeframe.  It  has  been  previously  stated  that 
the  routing  of  flows  for  the  entire  basin  is  accomplished  during 
each  timeframe.  During  each  time  interval  the  channel  flow  in  a 
upstream  reach  is  completely  combined  with  the  intervening  local 
flows  in  the  next  downstream  reach.  Bank  storage  or  time-lagging 
effects  are  not  considered.  For  monthly  timeframes  there  is  no 
need  to  consider  these  effects.  For  small  timeframes,  say  weekly 
or  less,  a  routing  technique  to  correctly  position  the  flows  in 
time  is  needed.  The  program  utilizes  the  modified  Gilcrist  flood 
routing  technique  to  accommodate  varyin°  times  of  travel  and. 
hydraulic  characteristics.  At  present,  there  is  no  data  input 
for  this  feature.  The  program  is  coded  to  set  the  coefficients 
equal  to  unity.  Coefficients  equal  to  unity  imply  complete  trans¬ 
lation  of  flows  from  one  reach  to  the  next  during  the  same  time 
period.  This  feature  can  be  generalized  by  removing  these  state¬ 
ments  from  the  program  and  by  inputing  the  necessary  coefficients. 
However,  the  number  of  timeframes  that  can  be  considered  is  twelve. 
This  limitation  is  binding  because  of  the  way  the  program  tallies 
deficiencies,  flows,  storage  levels,  and  pollutant  concentrations. 
Thus  the  program  could  consider,  for  example,  a  three  day  routing 
period  where  each  time  interval  is  equal  to  three  hours. 

8.  Waste  sites.  The  import  of  water  from  another  basin  may 
be  represented  by  a  waste  site. 

9.  Test  sites.  A  t*st  site  may  be  specified  to  receive  explicit 
and  conditional  pipeline  withdrawals  from  the  same  reservoir. 
Furthermore,  the  same  test  site  may  be  specified  to  call  for  further 
augmentation  by  channel  releases  from  up  to  ten  upstream  reservoirs. 
To  divert  water  from  a  stream  where  there  is  no  reservoir  (pump- 
pipeline),  place  a  reservoir  having  unit  capacity  at  that  location. 
This  is  limited  to  explicit  pipeline  diversions. 

10.  Reservoirs.  Reservoirs  in  series  cannot  be  effectively 
treated  by  this  program.  If  the  purposo  of  a  study  is  to  consider 
the  storage  requirements  for  satisfying  only  instream  flow  require¬ 
ments,  then  second  and  third  tier  reservoirs  may  be  specified  to 
make  conditional  pipeline  diversions  down  to  the  test  site.  The 
reservoir  on  the  highest  tier  is  first  depleted.  If  a  demand  still 
exists  at  the  test  site, releases  are  then  made  from  the  reservoir 
on  the  next  tier. 


71 


11.  Additional  notes  for  card  input  -  Option  It. 

1.  Card  ;  ir  this  program  is  used  at  an  installation  where 

magnetic  tapes  are  used  for  storing  intermediate 
results  (scratch  devices),  then  the  number  of  runs 
should  be  limited  to  eight.  This  constraint  is 
necessary  because  t lie  intermediate  system  response 
output  for  eight  runs  generally  fills  a  reel  of 
tape.  This  limitation  is  not  as  binding  if  scratch 
disk  devices  can  be  used. 

2.  Card  S-15  Often  it  is  desirable  to  have  the  historical  flows 

on  binary  tape  when  the  historical  trace  is  used 
to  test  alternative  plans.  This  saves  the  duplica¬ 
tion  of  decks  of  historical  flows.  A  binary  tape 
g  H-i" ..•‘I  f;  •  con  be  used;  however, 
this  tape  must  be  prepared  external  to  this  option. 

3.  Card  8-20  This  option  is  used  to  position  the  binary  flow 

tape  at  the  beginning  of  each  run.  It  may  be  used 
to  test  the  hydrologic  sensitivity  of  a  plan  by 
using  different  traces  of  operational  hydrology. 
Furthermore,  the  same  trace  can  ge  positioned  to 
test  different  plans. 


- 
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Output  for  Option  D 
Sample  Problem 


proposed  reservoir 


control  point  -  set  of  seasonal  flow  requirements 


gaging  station 


Given:  1)  Seasonal  coefficients  for  translating  flows  at  gaging 

station  to  inflows  at  dam  site  starting  with  the  month  of  April: 
0.292,  0.301,  0.235,  0.214,  0.279,  0.2S6,  0.299,  0.340,  0.3:9, 
0.30S,  0.272,  0.281; 

2)  Seasonal  coefficients  for  calculating  the  unregulated 
inflows  (between  dam  site  and  use  point)  from  the  inflows  at  the 
dam  site:  1.902,  1.835,  2.622,  2.995,  2.055,  1.855,  1.511,  1.519, 
1.812,  2.135,  2.042; 

3)  60  years  of  historical  flows  at  the  gaging  station; 

4)  Twelve  seasonal  flow  requirements  of  760  cfs-mo.  at  the 
use  point; 

5)  A  reservoir  whose  volume  is  equal  to  100  cfs-mo.  and  whose 
starting  contents  are  assumed  full  at  the  start  of  the  routing 
period;  and  that; 

6)  The  mode  of  reservoir  operation  is  to  always  meet  down¬ 
stream  flow  requirments. 

Find:  The  statistical  characteristics  of  deficiencies  at  the  test 
site  and  the  seasonal  probability  distributions  of  storage 
levels  in  reservoir. 
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Option  P  Output 

P»ge  1  Data  Input  echo 

a)  Seasonal  flows  at  the  dam  site. 

b)  Seasonal  coefficients  for  translating  flows  at  the  gaging 
station  to  inflows  at  the  dan  site. 

c)  Seasonal  coefficients  used  to  calculate  the  unregulated 
inflows  (between  the  dan  site  and  use  point)  from  the 
inflows  at  dam  site. 

d)  If  reservoir  operating  modes  2  or  3  have  been  specified, 
then  the  coefficients  for  the  set  of  regression  equations 
for  computing  seasonal  releases  will  appear. 

Page  2  Statistics  of  historical  flows  at  dam  site. 

a)  Seasonal  means,  standard  deviations,  skewness  coerficients, 
lag  one  correlation  coefficients,  and  the  coefficients 

of  variation. 

b)  Grand  seasonal  mean,  standard  deviation,  lag  one 
correlation  coefficient,  and  coefficient  of  variation. 

c)  Annual  mean,  standard  deviation,  skewness,  lag  one 
correlation  coefficient,  and  coefficient  of  variation. 

d)  Seasonal  histograms  of  flows. 

•)  Histograms  for  all  seasons. 

f)  Histograms  of  annual  flows. 

The  output  for  the  histogram  analysis  is  intrepreted  in  the 
folllowing  manner:  Consider  October:  1)  7/60  or  11.6%  of 
all  Octobers  the  flow  is  equal  or  less  than  SO. 28  cfs; 

2)  1S/60  or  25%  of  all  Octobers  the  flow  is  equal  to  or 
less  than  7S.42  cfs;  and  3)  48/60  or  60%  of  Octobers  the 
flow  is  equal  to  or  less  than  S02.8  cfs. 
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Three  different  transformations  of  the  flow  data  are  tried  with 
the  purpose  of  finding  the  transformation  which  minimizes 
the  skewness  of  the  data.  The  three  transforms  are:  no  trans¬ 
form,  log  10,  and  square  root.  An  overall  skewness  coefficient 
is  determined  by  weighting  the  seasonal  skewness  coefficients 
by  the  means  of  their  respective  seasons.  The  three  diagnostics 
that  appear  are  the  weighted  coefficients  for  each  transform 
divided  by  the  standard  deviation  of  the  skewness  coefficient. 
This  standard  deviation  is  computed  from  a  formulation  by 
Kendall  and  depends  only  on  the  sample  size  (5).  The  trans¬ 
form  whose  diagnostic  is  minimal  is  chosen  to  be  used  later 
in  the  streamflow  simulation  analysis. 

Pages  4-11  are  devoted  to  describing  the  statistical  performance 
of  the  system  using  the  historical  period  of  record.  Pages 
11-19  are  the  exact  counterparts  using  a  trace  of  operational 
or  synthetic  hydrology.  They  will  not  appear  nor  be  discussed. 
However,  the  print  suppress  features  for  the  synthetic  record 
analysis  output  will  be  indicated  in  the  discussion  of  the 
historical  record  routing  output. 

Page  4 -  Statistical  moments  of  deficiencies. 

a)  Shown  for  each  season  is  the  probability  of  failure,  the 
mean  of  all  the  failures,  the  standard  deviation,  the 
maximum  miss,  and  the  minimum  miss.  These  failures  are 
tallied  without  regard  for  the  duration  of  the  miss. 

Section  c  deals  with  these  statistics.  On  input  card  7 
for  this  option  a  non-zero  punch  in  column  2  will  suppress 
this  output  in  the  historical  record  analysis.  A  non-zero 
punch  in  column  14  will  suppress  this  output  in  the 
synthetic  record  analysis. 

b)  The  annual  probability  of  failure,  the  mean  annual  failure, 
the  standard  deviation,  the  maximum  annual  failure  and  the 
minimum  annual  failure  ore  then  presented.  An  annual 
failure  is  tallied  if  the  flow  requirement  for  any  season 
during  that  year  is  not  maintained.  The  magnitude  of  an 
annual  miss  is  the  sum  of  all  seasonal  misses  during  that 
yaar.  A  non-zero  punch  in  columns  4  and  16  will  suppress 
this  output  in  the  historical  record  and  in  the  synthetic 
record  analysis  respectively. 
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c)  The  probability  of  failure,  the  mean,  and  the  standard 
deviation  for  failures  of  different  durations  are  shown 
A  duration  of  a  continual  failure  may  be  as  long  as  36 
seasons.  The  magnitude  of  a  miss  for  any  stated  duration 
includes  all  prior  consecutive  misses.  A  non-zero  punch 
in  columns  6  and  18  will  suppress  this  output  in  the 
historical  and  in  synthetic  record  analysis, respectively. 

Page  5  Distribution  of  failures. 

The  accumulative  histograms  of  failures  for  each  season 
and  for  the  annual  failures  appear  in  tabular  form. 

Consider  the  month  of  September;  hal^  of  the  misses  were 
equal  to  542  cfs  or  less.  Consider  the  annual  misses; 

95%  of  the  annual  misses  were  equal  to  3283  cfs  or  less. 

A  one  punched  in  columns  8  and  20  will  suppress  this  output 
in  the  historical  record  and  in  the  synthetic  record 
analysis  respectively. 

Page  6  Distribution  of  storage  levels. 

The  probability  distributions  of  storage  volumes  for  each 
season  is  shown  in  tabular  form.  Consider  the  month  of 
Juno;  in  33%  of  the  Junes,  the  storage  in  the  reservoir 
was  equal  to  10  cfs-mo.  or  less.  For  the  month  of  March, 
the  storage  in  the  reservoir  was  always  between  90  cfs-mo. 
or  full.  A  one  punched  in  columns  10  and  22  will  suppress 
this  output  in  the  historical  record  and  in  the  synethetic 
record  analysis  respectively. 

Moments  of  flows  past  c,est  points. 

The  mean  and  the  standard  deviation  for  the  flows  at  ihe 
test  site  are  shown  for  each  season.  A  non-zero  punch  in 
columns  12  and  24  will  suppress  this  output  in  the  historical 
record  and  in  the  synthetic  record  analysis  respectively. 


Pages  7,  8,  and  9  Accumulative  histograms 
failures. 


>f  seasonal  and  annual 


Presented  on  page  7  in  graphical  form  are  the  accumulative 
distribution  of  failures  for  the  first  six  months  of  the 
year  used  in  the  analysis.  Here  the  first  month  is  April. 

The  remaining  6  months  will  be  found  on  page  8.  The  ordinate 
scale  of  the  plot  is  the  accumulative  probability;  each  tic 
being  equal  to  .02  or  2  percent,  (hi  the  absicca  the  scale 
is  the  percentage  of  the  range  or  maximum  miss  expressed  in 
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decimal  form.  The  maximum  failures  for  each  of  the  six 
seasons  are  shown  directly  below.  The  numeric  entries 
(••l”  -  "6”)  are  plotted  at  10  percent  intervals.  In  each 
block  they  appear  spaced  from  left  to  right  ("6"  -  "1")* 
always  to  the  left  of  their  grid  point.  This  constraint 
is  necessary  to  preclude  the  possibility  of  overlap. 

Simply  displace  each  entry  to  the  nearest  10th  point  to 
the  right  and  read.  Consider  the  following  examples. 

For  the  month  of  July  and  August,  6  percent  of  the 
deficiencies  were  equal  to  6S  cfs  or  less.  For  the  month 
of  June,  80  percent  of  the  misses  were  equal  to  (.7MS6.) 
or  319.2  cfs  or  less.  Shown  on  page  9  is  the  accumulative 
histogram  of  annual  failures.  A  two  punched  in  columns 
8  and  20  will  suppress  both  the  tabular  and  graphical 
displays  in  the  historical  record  and  in  the  synthetic 
record  analysis  respectively. 

Pages  10  and  11  Storage  level  probability  distributions. 

The  distributions  of  storage  volumes  appear  in  graphical 
fora  for  the  first  six  months  of  the  year  on  page  10. 

The  distributions  for  the  remaining  6  months  appear 
on  the  following  page.  The  ordinate  scale  of  the  plot 
is  the  probability  of  occurrence  and  not  the  accumulative 
probability.  The  absicca  scales  the  reservoir's  capacity. 
Consider  the  month  of  June.  The  storage  in  the  reservoir 
was  10  cfs-mo,  or  less  for  34  percent  of  the  Junes.  A 
two  punched  in  columns  10  and  22  will  suppress  the  tabular 
and  the  graphical  displays  in  historical  record  and  in  the 
synthetic  record  analysis  respectively. 

Page  20.  Statistics  of  generated  flows.  The  statistical 
parameters  for  the  first  100  years  of  the  generated  trace 

.  at  the  dam  site  appears  on  this  page.  If  less  than  100  years 
were  specified,  then  the  entire  trace  will  be  used  to 
compute  the  parameters.  This  output  is  precisely  like 
that  presented  for  the  historical  trace  on  page  2.  A  non¬ 
zero  punch  in  column  26  will  suppress  this  output. 


W  KtStKVtlU  CAPACt  11  capacity  *  too. 
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Output  for  Option  NMP 
A.  Historical  trace 

1.  Seasonal  ir^ans,  standard  deviations,  and  skewness  coefficients 
for  each  site, 

2.  Grand  means,  standard  deviations,  and  skewness  coefficients 
for  each  site. 

3.  Skewness  transform  diagnostic. 

4.  Spatial  correlation  matrix  (M0). 

5.  Time-space  correlation  matrix  (M_j).  The  diagonal  elements 
in  this  matrix  are  the  lag  one  correlation  coefficients 
for  each  site.  The  off-diagonal  elements  in  this  matrix 
are  the  lag  one  correlation  coefficients  of  each  site 

with  all  of  its  neighbors. 


6. 

Intermediate  matrices:  A, 

Simulated  trace 

(1st  100  years). 

1. 

As  above. 

2. 

As  above. 

3. 

Matrices: 

a)  M 

LX  1IU 

b)  M_  j 
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Reservoir  capacities  arc  200,  100,  300,  and  400  at 
coordinates  1,  3,  8  and  11, respectively. 

The  starting  contents  of  the  reservoirs  are  100,  50,  100, 
and  200,  respectively. 

There  are  two  waste  sites  in  the  system  located  at 
coordinates  6  and  9.  The  pollutant  flow  rate  at  the 
first  site  six  is  50  and  its  concentration  is  1200. 

The  pollutant  flow  rate  at  coordinate  D  is  100  and  its 
concentration  is  300. 

The  test  sites  at  coordinates  4,  2,  5,  and  13  have 
flow  requirements  of  1000,  250,  350,  and  1000,  respec¬ 
tively.  The  test  site  at  coordinate  7  has  a  water 
quality  standard  of  75. 

The  reservoir  at  coordinate  1  may  make  channel  releases 
to  maintain  the  flow  requirements  at  the  test  sites 
located  at  coordinates  2  and  5.  The  storage  in  the 
reservoir  at  coordinate  2  is  to  be  kept  half  full. 

The  reservoir  at  coordinate  8  may  make  channel 
releases  to  maintain  the  water  quality  standard  at 
coordinate  7.  The  reservoir  at  coordinate  11  may  make 
8  pipeline  withdrawal  to  maintain  the  flow  requirement 
at  coordinate  13. 

The  concentration  of  the  initial  storage  in  all  reservoirs 
is  equal  to  zero. 

The  unregulated  inflow  for  coordinates  1-5,  6-9,  and  10-14, 
arc  defined  by  scaling  the  flows  at  gaging  stations  1,  2, 
and  2,  respectively.  The  scaling  factors  for  the  U 
coordinates  arc:  1 . ,1. ,1 . ,0. ,1 . ,1. .0. ,1 , ,1 . ,i . ,1 . ,1. , 

1.,  and  0.  These  factors  multiplied  times  the  gage  flow 
define  the  net  inflow  into  each  coordinate.  This 
corresponds  to  choice  A  in  card  *6. 

The  background  concentration  of  all  the  unregulated  inflows 
is  equal  to  :cro. 


Find:  For  a  ■'•outing  period  of  one  season,  determine  the  concentration 
at  coordinate  7,  the  storage  levels  in  all  reservoirs,  and  the 
tally  o:  deficiencies  at  all  test  sites.  The  gaging  station 
flows  j.re  100,  200,  and  300,  respectively. 

The  next  two  pages  arc  a  listing  of  the  card  input  necessary 
to  study  this  problem.  The  ordering  of  the  test  and  junction 
points  should  be-  carefully  noted.  The  test  point,  at  coordinate  4 
is  the  first  test  site  in  the  system.  The  test  point  at  coordinate 
13  is  the  5th  test  site. 


OUTPUT 

Page  1  Data  echo.  Includes:  a)  title  of  study;  b)  years  of  routing; 
c)  number  of  pollutants;  d)  number  of  years  the  flow  tape  is 
to  be  skipped  forward;  e)  random  number  generator  trigger; 
f)  starting  point  in  system;  g)  choice  for  computing  unregulated 
inflows;  and  h)  tabular  map  of  system. 

Pag  2  Data  echo.  Description  of  each  coordinate. 

Page  3  Data  echo.  Includes:  a)  coordinate  to  gaging  station 
reference  map;  b)  gaging  station  scaling  factors  for  each 
coordinate;  c)  routing  coefficients  (set  equal  to  one  in  the 
program);  d)  reservoir  capacities,  initial  storages,  and  dead 
storages;  and  e)  parameters  of  the  seasonal  operating  rules 
for  the  reservoirs. 

Page  4  Data  echo.  Includes:  a)  number  and  location  of  reservoirs 
that  may  make  channel  releases  to  flow  requirements  or  water 
quality  standards  at  specified  test  sites;  b)  allocation 
option  -  node  of  calculating  releases  for  more  than  one 
reservoir  serving  a  test  site;  c)  test  site  (number  and 
location)  that  may  receive  pipeline  supplements  from  dam 
site  (number  and  location);  d)  evaporation  of  storage  in 
reservoir  option;  and  e)  concentration  of  starting  storages. 

Page  5  Data  echo.  Includes:  a)  minimum  waste  flow  rate  and  minimum 
poTlutant  weight  for  each  waste  site;  b)  seasonal  means  flow 
rate,  standard  deviation  of  flow  rate,  seasonal  mean  waste 
weight  for  each  pollutant,  seasonal  standard  deviation  of 
waste  weight  for  each  pollutant  (note  -  for  the  1st  waste  site 
at  coordinate  6,  the  waste  weight  is  30000(100*500));  seasonal 
Kean  and  standard  deviation  of  flow  requirements  for  each  test 
aitc;  *nd  d)  seasonal  water  quality  standards  for  each  test  site. 
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Pr,ge  6  Tally  of  deficiencies^  On  the  left  half  of  the  page  is  the 
number  of  seasonal  deficiencies  for  each  test  site  in  the 
system.  A  deficiency  is  tallied  whenever  any  of  the  water 
water  quality  standards  or  flow  requirement  at  a  test  site 
has  not  been  maintained.  On  the  right  half  of  the  page  is 
the  number  of  times  that,  the  critical  or  greatest  failure 
was  attributed  to  the  water  quality  standards  or  flow 
requirement.  The  count  is  independent  of  season.  In  the 
second  line  there  is  an  entry  entitled  "shortage".  This 
is  a  count  of  the  number  of  times  when  there  was  a  deficiency 
at  a  test  site  and  that  there  was  available  storage  in  any 
of  the  specified  reservoirs  to  satisfy  this  deficiency. 

This  occurs  when  reservoirs  in  scries  have  been  specified 
to  satisfy  by  channel  releases  downstream  test  sites. 

Page  7  Probability  of  misses.  The  probabilities  of  seasonal 

HefTcicncies,  critical  failures,  and  yearly  failures  are 
presented.  A  yearly  failure  is  counted  when  the  water 
quality  standards  or  flow  requirement  for  any  season  during 
a  year  is  not  maintained. 

If  all  the  print  suppress  options  have  been  depressed,  then 
the  output  terminates  with  page  7.  The  output  for  the 
second  run  in  a  multi-run  job  begins  with  page  2;  that  is, 
all  the  design  information  is  reprinted. 

Page  8  Seasonal  deficiency  summary.  The  number  of  misses,  the  mean 
of  the  misses,  the  standard  deviation,  and  the  probability 
of  failure  is  presented  for  each  water  quality  standard  and 
flow  requirement.  A  water  quality  standard  deficiency  or 
miss  is  the  amount  of  dilution  water  needed  to  lower  the 
test  site  concentration  down  to  the  level  of  the  standard. 

The  calculation  assures  that  the  concentration  of  the  diluting 
water  is  equal  to  zero.  Later  output  will  give  the  means 
and  standard  deviations  of  the  pollutant  concentrations  in 
the  reservoirs.  Non-zero  punches  in  columns  40,  45,  and  50 
in  the  first  card  of  *9  and  in  columns  5,  10,  and  15  of  the 
second  card  of  *9  will  suppress  this  output. 

Page  9  Maximum  seasonal  deficiency  summary.  This  is  a  summary  of  all 

the  maximum  seasonal  failures.  The  number  of  failures,  the  mean, 
the  standard  deviation,  the  probability,  and  the  maxima  is 
presented.  Note  that  the  sum  of  all  the  seasonal  failures 
shown  on  page  6  for  a  test  site  must  equal  the  number  of 
deficiencies  shown  here,  A  non-zero  punch  in  column  20  in 
the  second  card  of  *9  will  suppress  this  output. 
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Page  10 


Page  11 


Yearly  deficiency  summary.  The  number  of  yearly  failures, 
the  mean,  the  standard  deviation,  the  probability,  and  the 
maxima  of  the  yearly  failures  is  shown  for  each  test  site. 
A  yearly  failure  is  the  sum  of  all  the  maximum  seasonal 
failures  during  that  year.  A  non-zero  punch  in  column  25 
in  the  second  card  of  *9  will  suppress  this  output. 

Reservoir  performance  summary.  Show  on  this  page  are  the 
numoer  of  times  the  reservoirs  spilled  and  the  number  of 
times  they  were  depleted  down  to  dead  storage  or  emptied. 

A  spill  is  tallied  when  the  total  assets  (inflow  plus 
start  of  the  season  storage)  exceed  the  sum  of  all  the 
withdrawals  and  releases  and  the  capacity  of  the  reservoir. 
A  non-zero  punch  in  column  30  of  the  second  card  of  *9 
will  suppress  this  output. 


Page  12  Reservoir  performance  summary  -  distribution  of  storage  levels. 
The  mean  and  the  histogram  of  storage  volumes  for  each 
reservoir  for  each  season  is  presented.  The  capacity  of  the 
reservoir  is  the  range  of  the  histogram.  If  twelve  seasons 
had  been  specified,  then  11  more  pages  would  follow.  A  non¬ 
zero  punch  in  column  30  in  the  second  card  of  *9  will 
suppress  this  output. 

Page  13  Distributions  of  maximum  seasonal  deficiencies.  The  maximum 
seasonal  failure,  the  accumulative  histogram,  the  mean,  and 
the  standard  deviation  of  the  maximum  failures  for  each 
season  is  presented  for  each  test  site. 


The  range  of  the  histogram  is  the  maxima  seasonal  failure.  If  12 
seasons  had  been  specified,  then  11  more  pages  would  follow. 

A  non-zero  punch  in  column  35  in  the  first  card  of  *9  will 
suppress  this  output. 

Page  14  Distributions  of  maximum  yearly  deficiencies.  The  accumulative 
Histograms  "of  annual  fai lures  are  shown  for  each  site.  The 
range  of  the  histogram  is  the  maxima  of  the  yearly  failures. 

A  non-zero  punch  in  column  25  in  the  second  card  of  *9  will 
suppress  this  output. 


u»  mm&l 


Page  15  Histograms  of  deficiencies  for  water  quality  standards.  The 
histograms  of  deficiencies  for  water  quality  standard  1  are 
shewn.  The  range  of  each  histogram  appears  Jirectly  below, 
(f  12  seasons  had  been  specified,  then  11  more  pages  would 
follow.  If  more  than  one  pollutant  had  been  specified,  then 
the  series  of  pages  would  be  repeated.  Non-zero  punches 
in  columns  5,  10,  15,  and  20  in  the  first  card  of  *9  will 
suppress  this  output. 

Page  In  Histograms  of  seasonal  deficiencies  for  flow  requirement. 

The  histograms  of  deficiencies  for  the  flow  targets  at  each 
of  the  test  sites  arc  shown  on  this  page.  The  range  of  each 
histogram  is  the  maximum  seasonal  miss  and  is  shown  directly 
below.  If  12  seasons  had  been  specified,  then  11  more  pages 
would  follow,  A  non-zero  punch  in  column  30  of  the  first 
caid  of  *9  will  suppress  this  output. 

Pave  17  Seasonal  statistics  of  system  flows.  The  maximum  flew,  the 
Tu s tog rams ,  the  mean,  and  the  st; noard  deviation  of 
the  seasonal  flows  for  each  coordinate  are  presented.  The 
range  of  the  histograms  are  the  maximum  seasonal  flews.  At 
test  sites  the  lower  fifth  of  the  histogram,  that  is,  the 
flows  between  zero  and  0.2  of  the  seasonal  maximum,  are 
further  subdivided  into  fifths.  This  is  designated  by 
"LON  FIF."  at  each  of  t.ie  test  and  junction  points.  If  12 
seasons  had  been  specified,  then  11  more  pages  would  follow. 
A  non-zero  punch  in  column  35  in  the  second  card  of  *9  will 
suppress  this  output. 

Grand  means  and  standard  deviations.  The  mean  and  standard 
devi'  an  for  all  seasons  is  presented  for  each  test  site. 

A  rev  ro  punch  in  column  3.5  in  the  second  card  of  *9  will 
suppress  this  output. 

System  concentration  summary.  For  each  test  site  the  mean 
amT  the  star  lard  deviation  of  the  system  concentration  for 
the  pollutant  is  sno.vn.  1/  12  seasons  had  been  specified, 
then  11  more  columns  of  means  and  standard  deviations  would 
appear.  If  more  than  one  pollutant  had  been  specified,  then 
another  page(s)  would  follow.  A  non- zero  punch  in  column  40 
in  the  second  card  of  *9  will  suppress  this  output. 

Page  20  Distribution  of  system  concentrations.  The  minimum,  the 
FTst og ram,  'arid  “the  maximum  Vf  the  system  concentrations  is 
presented  for  each  test  site.  The  range  of  the  histogram 
is  the  maximum  concentration  less  the  minim1  .*  concentration. 
The  histogram  is  formed  over  all  seasons.  If  more  than  one 


Page  IvS 


Page  19 
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pollutant  were  specified,  then  another  pagc(s)  wc  ^ld  follow. 

A  non-zero  punch  in  colum  40  in  the  second  card  of  *9  will 
suppress  this  output. 

Page  21  Concentration  of  pollutants  at  reservoirs.  Shown  for  each  res¬ 
ervoir  is  the  mean  and  the  standard  deviation  of  the  pollutant, 
for  all  seasons.  If  more  than  one  pollutant  had  been  specified, 
then  the  mean  and  standards  deviation  for  each  site  would  be 
grouped  together.  A  non-zero  punch  in  column  40  in  second 
card  of  *9  will  suppress  this  output. 

Page  22  Distribution  of  concentrations  at  the  dam  sites.  The  minimum  con¬ 
centration,  the  histogram,  and  the  maximum  concentration  is 
siiown  for  each  reservoir  site.  The  histogram  is  formed  over 
all  seasons  and  its  range  is  the  maximum  concentration  less 
the  minimum  concentration.  Histograms  for  other  pollutants 
would  appear  on  other  pages.  A  non-zero  punch  in  column  45 
in  the  second  card  of  *9  will  suppress  this  output. 


I 

t 

j 

i 

i 


APPENDIX 


Listing  of  Option  D. 
Listing  of  Option  N'MP. 
Listing  of  Option  R. 


DEM 


. . .  (Main) 
Subroutines 


\ 


READ 

STAT 

SEAMON 

PRINl 

XBARSD 

XSKEW 

XCOR 

HIST 

TRANS 

SETUP 

CLT24 

GBNFLO 

BACTAN 

ROUTE 

PROCES 

FAILDR 

ABTRAC 

OUTPUT 

OUTPIT 

FIRST 

PACE 


PRTLOT 


HUH  f.t-oiPBATUW  CUHTKfH-S.  AM>  flt«T  1 1 IWALL  V  .  STa 


PKOGKAN  CVCtE  POINT  FUR  NEW  SET  OF  HAMS  AND  niLOTIIIN  REUOS 


•  • 

« v  *». 


tITlltG.il 


.Y:- 


Ml 


si 


f; 


•  * 
►  • 


i/I  • 

■—  til 

S  *s 


s  § 


e  • 

Ml  X 
►*  ** 

«  9 

*  O' 

M.  • 

«s  o*4 
%  «1 
3».  3  *3 

g?  2  *1 

*«•*»>! 
Off  <1  ►  «4 

6  —  *-  j  ■*  ►* 

-*  —  w»  ii»  5 

-  x  1 «  a 

3  — 

5X  5 

*—  x  ,  -*  ** 

•  —  I  i  -*  *• 

St  a  IS5 


<r  c  ~ 

fNi  fO 
(fl  s*\ 


*,0*  C1-*  Oilm  *  ICC 
m  ft  #r. 


11 


« 

2 

3  vw 

5  ! 

*  i 


uu  >M« 


STAT(FLUHtS,*VRS,N*€AS.  IO  1ST  « 1  STRUT  ( I  FA W  ,  XM  •  Sl)t  COR  t  REfi  451 


SWAM*  l-A!  Ol  INI 


••.X-M.s  M, 

<< 


■.  •  »  V  V-TS 

V'.<  A? 


kmWWMMII  - 


■«  V;pv.V^ 

*  r*  ; 


* 

'£rf- 

9»- 
"1 


■fimi 


in 


r«n  m  ^  ir  <  r^aqO'd-'i 
oc  oqcc  c  q  o-4«  , 
3  3 


«-«  pJ  m  *1  ir\  <  rp-  oc  3>  c  <-<  no)  m  >r|  x'  <c 

<*>  ifj  m  tri  &.  in  tin  rr,  w  <  •»■ 

<0  <*  -Cl-C  C  *  -Cj-C  <  *  .£ 


•  « 
—  Z 
Q  « 


! 


o  j 

jt !«•  £ 

■«  2 

•  ,► 

2  3  i 

uu  */»  u 


'7 

5-5x  = 

ac  *  ui  r  x 
4L  ►» 

C‘  -I  *  - 
«u<5i 

UA,UW,W 

N  * 

ivj  iNi  t 

i 


► 

* 

•z' 

«  * 

in  (/« 

¥*  ,  JM 


x  !  c 

✓>  i  * 


&  > 

a.  v^r  -k  •  *» 

-  : 

s/l  X  3  ,'JJ 

42  XKs*  .VU 

<  *g>-  (a: 

♦  <  -«*  r 

-*  ft  — »  ►- 

w  j  ^  ; 

—  i  ’  t/i 

SJI  £ 

c?  =f^  « 


-  5  Ji  ■  ’  i 
!S  !i5S  ° 
*  ,-  :5  i'r  | 

W  it  '  <rt  J1.V  w 


l(0.lifti'H[S(IY.lSF*S» 


CnMPtiT  LUC 


753 


-L»mtH»i&t&.PAJICtMt»CKt*S£A«CRliyR»LHIST  .  INOFX  .  f)AT  I  .FMAX.FMA 
2FHAXVK , (HIST ) 

_  MlHfcaSfU*  PA-.XaUi.il  . P AR  S11U2.1U EA Ri II 12 . 3 ) .LH1 ST t  U  . ?n > 


D  17 


*-»o4  m  *  • 

OQOdO 

ttictqcci 


crc'*-«cJm«^,.r*cjr**x'j'C  <\j  <*s  tr  c  r*  ®  c!  **  fs<  •'*'  4I  *d  p* 

*«  oc  *l«  *jx  X  x  a  x  *,*  «.«  x  *  * k  xj*  X  at  *k  ««  x  «  «j«  x  *  «  *  *  «  *k 


»  p*  ad  »  Cl 

1  ■«•  •#  >»•  it; 

tdxi 


< 

«*  i 

«•  h 

5  * 

\  % 
—  < 
1* 

•  < 
M  k 

r  2 

S  j 

3  2 


+  i* 
O  ^ 


3  S  IS 
5  ^  L* 


*  1  <  - 


N  •  «/> 

«*  V*  — 
M  *-.  *  , 


•  *  1  >«l»l 

-  -  <  J  JEN  ■ 

-  O’  i  —  O 

-  -  x  3  -1*  . 

*•  >  —  ,  u»  *  — 


—  ,\M*« 

—  uJ  2  «- 
>■  U  l  “1- 
<  <  —  u  4 


!;  *  *1 

.  X  « 

I  5  'A 

j  #.  «r  » 

MMkv) 

<  —  I*  • 

ViJ  —  — 

s;:t 

•  <  vs 

«■*  UJ  *  4 

tt>  VS  £  W 

t/l  m  »  4* 

«  f*  M 

Ul  3  *  — 
—  vs  vs  x  * 

1 

a  >»  4  4»  w 


$  £  rgSES’-ST-"*-* 

I  |  ^  < »  H  =  *  3  *;  *  2  >  *  x 


:  ’  s  x  : 

i  x  x  O  *  a. 


.  *  vi 

a-  ’«* 

Jx|3  * 

«u«i 

1*4-3 

33153 

1-4  VS  X  ■  •-»  tm 

S  I  7  -  r  - 

§Si*Sf 

x|  X  U  X  X  * 


53S 

•  Of  f>4 
*  >  * 

>  XI  p* 
CMI 
at  ^ 

O  4  • 

*  Ml  -* 

at  ©  1. 

►  *■* 
tt  •»  * 
VS  *  j  iJ 


VS  UV  * 

•  •  VS 

at  «  at 
►  v*  ► 
<  C  2 

- 1  - 

tT  — ,  « 


c  =  < 

l<2  *  Z 

i^juj  x  7i 

3«/>  <  < 
X  UJ 
X  >A 


Xll 

<  I-  — 
X  C  X 
-  « 
UL  x  X 


C  «rt 
O  C  u. 


J  •"  z 

t/t  i/I  UJ 

X  —  *  o 

q  z  a 

-J  i/I  X  X  * 

£1  z  ?al 

J*  =  5” 
i  =  c 

CJojx: 

v*  :c 

5iu.ni. 

UJ  *  ■ 
vth-  c  Z 

z!: 

uJ  —  vs  • 

C  V>  c  < 

^  <  a-  i 

v>  a  v) ; 

z$"l- 

atJu.  S 

*  S' 

hZ  J  \, 

v?  —  I  x  < 

"•  o  ^  s 

X  <  *  Z  ' 
r  a  _i 
N/  1  li 
(T>  — 

©  X  >  J 
X-  «  < 
-<  x  *  r  • 
-  ©  c  i  ■ 


o  •  • 

—  • 

-X  < 
x  —  ( 

4-  ■!*> 

—  X  >  ■ 

C  X  •-  x 

iv  •  <  ; 

•>.«*■ 

—  >.—  u. ' 

3  <y.  c*  »  : 

X  I-—-I' 

•  z  -< 

«XXXI 

X  >  *  <  i 

•  (II  |  X  1 

x  x ' 


H  -H-  r  H*  * < 
:  ij<: 

il-  X  x.  X  <  *  4.  -j  1 


r  jj  «  m 
•S  —  X  fc. 
«  -  *  - 
>  J.  » 

Jill 


■  x  «  <  1/1 ; 

.  J  it  s;  z  : 

■  —  vs  Z:  •  i 

•  •  Z*-*« 

:  —  j  <  »  v 

•  Z  u  vs  < 

►  *  <’L'< 

•  -»vt : 

:  —  ac  vs  < 

•  1:1- 

>  *  UL  <  •  S 

t  •  *■  — 

•  -  X  a  vs  : 

*¥  <  w  4  : 
S  .  *  >-  UJ 

•  m  •  j  vs  vs 

•  w  1  ^  1 

I K ( 

•  VS-C4' 
l*  XJf  Z  ■ 

z::i 

-  C  "  C  1 

•  H  at  M  *  1 

■  <«J  < 

■  —  r  *  —  1 
I  1  4  ii  4 

i9i~- 

:  »  ;  k  * 

.  X  U.  -  X  1 


XU.  4 
X 

1*1  UJ  * 

I*  c  c 
c  z  **■ 

X  « 

-4  X  'X 
-  X  c 
U’«tl 
<4  4  X 


3|x82^ 


u  uuu 


1ftF-JL.0i2j.2uJ. 


soev-o 


X*TAK-  seasonal 


UMftEftULATen  FLOWS  <H6TWeeW 


I mast  t - 1 

X"M  »*I  ISRASJ— *«UKU  SI4\I_ 


MEANS 


•  *'<«£- 


l«VSAV  ♦«lt**XX 


starting  storage  of  reservoir 


ptukn  fuims  h flow  hfsfrvoik 


CONTINUITY  check 


no*-  nows 


AS1 


SUHMUUTIME  *AlLMUHMMfXA!SS, 
CUMMir  IS  ccwsecuf  IVE  .fAltUftfc- 


1 59  A 


St  ASUNAL  AMD  ANNUAL  HISTOGRAM  PkOMAH IL  HIES 

. . . . . . IMh. 


Hi 


oouU  uyO'joaa 


!  I 


|  : 

,L 

r-  a 

C*  O 

m  < 

tr  < 

P-  X 

O'  c 

(NJ 

IT  ^ 

!*•  3C 

O'  c 

-« 

m 

cr>  ^ 

f“*  CQ 

*c 

^  C\i 

** 

If  »£ 

f-  X 

O'  O 

-N 

PO  4- 

tf*  <\  ! 

i 

{ 

9  0 

O'  a  c  c 

c  cf 

c  o  c  c 

c  -4 

•— »  4-4 

•— «  *- 

1— X 

*— <  \ 

fV  p\ 

(NJ  (\J 

f\ 

PSi  rv 

(X  rr 

(T 

m  m 

PO  <r. 

m  fr. 

ro  >t 

^  4T 

■4* 

-4  -r 

*  • 

« -t 

^  p-* 

p-  r- 

r-  p> 

p-  p- 

p-  Pi 

r-  p- 

r-  c- 

r- 

p-  r- 

r- 

J^4 

-  K 

r-  r^t 

r-  p- 

t*“  p- 

r- 

r-  Pi 

p*  f- 

P-  p* 

P“  ^ 

p-  p4 

P“  P-1 

! 

r 

#■ 

■■ 

H  ^4 

H  <— 

*—*  «*» 

4—4  4— 

#— t  »»  ’ 

4—4 

rH  4-4 

r-4  4—1 

■ 

«— 4  »— 

4H  4— 

—4  ♦— 

4*4  —-4 

•—4  -H 

^  4  4 

*-*  4—4 

4~ 

f— <  ^-4 

p 

c 

£ 

X 

r 

! 

X 

< 

VC 

1  i  ! 

. 

< 

o 

VC 

i 

i 

-J 

* 

X' 

1  : 

G 

X 

X  X 

• 

X 

• 

<  Q. 

>• 

_ 

< 

cr 

-j  a 

8 

< 

h- 

« 

D  3 

l; 

G 

4- 

a  VC 

1 

2 

• 

< 

» * 

7 

CO 

X 

p-  X 

J 

-J 

> 

X 

l  1, 

CL 

• 

< 

VC 

t  '* 

•f 

> 

VC 

V" 

'JJ  — 

f  1 

CL 

X 

VC 

i 

«/ 

-J 

G 

VI 

G 

f  i 

i 

_ 

£ 

•J 

Z 

<  V^ 

t  1 

1 

— 

<1 

X 

X 

X 

O  u- 

►- 

L 

od 

¥-4 

p“ 

< 

: 

z 

c 

C 

c 

I  3 

V 

.  1 

►- 

c 

> 

c~ 

•x 

W  ^ 

X 

■  ! 

* 

a 

H" 

X 

z 

p- 

X 

5 

VC 

a 

3 

V 

UJ 

< 

— 

3  i 

Vi 

1 

i  ' 

H 

•H 

CO 

G 

X 

< 

X  P* 

•J 

►- 

2 

i 

CL 

a 

X 

1  3 

1 

■ 

1 

K* 

vt 

> 

X 

• 

z 

Cf 

—  X 

! 

• 

Z 

CL 

V 

V 

u;  -C 

x- 

_ 

«J 

.  1 

»- 

V 

e. 

v.  * 

c 

X 

•  VC 

< 

3 

4 

• 

»— 1 

G 

G  • 

VC  X 

3 

'  :  ! 

a 

P- 

. 

£ 

Z 

CL 

c 

Z  V" 

z 

X 

VC 

< 

3 

VC 

— > 

2 

X  z 

a 

a 

V' 

< 

<- 

VL 

h 

Z 

vj  C 

•— 

c  < 

< 

i  ) 

VC 

A. 

X 

z 

< 

C  P- 

j- 

pu 

X  D? 

{ 

X 

i 

S 

VI 

CL 

X 

X  ►- 

■j 

Z 

q 

I 

V* 

s 

» 

'JJ 

> 

c 

< 

X 

G  X 

z. 

i 

p- 

V? 

v*  u. 

G 

vr  C 

X 

X 

<i 

X 

• 

vc 

J 

1/ 

G 

-1 

VI 

CL  | 

>— 

3: 

►4 

X 

p  z 

G 

■ 

< 

■H 

a 

CL 

M  V 

VC  — 

VC 

G 

X 

r  i 

CL 

CL 

X 

a 

X  CD 

VC  Z 

X 

— j 

X 

VC 

z  a 

i 

at 

• 

V 

4—4  G 

X 

J.  VI 

G 

vc 

—  1 

z 

*  , 

VC 

« 

> 

z 

Z  v' 

X 

■« 

X 

X 

r\ 

a 

X  . 

» 

Z 

x 

X 

> 

c  - 

u. 

-j 

CL  CO 

X 

c 

c  1 

5 

> 

* 

* 

UJ 

vi  X 

z  < 

< 

G  > 

VC 

OL 

X 

< 

O  1 

j 

mJ 

X 

• 

«> 

G 

< 

-4  Z 

X 

z 

•  ^ 

u 

X 

>- 

\U 

— •  1 

x 

c 

a. 

X 

x 

UJ  -J 

••  — 

< 

v*  < 

3  Z 

i 

X  ' 

{ 

o 

• 

i 

WM 

VC  « 

l»-  / 

•J 

VC 

4* 

V 

a  x 

1 

i 

r  \ 

Z 

» 

z 

p- 

z 

G 

» 

’JW 

•  ~ 

<  < 

3 

< 

.1  < 

VC 

P»  V* 

•ij 

« 

a 

u 

4 

it 

*| 

111  f 

X  v. 

Jl  e 

u. 

— 

-J 

3  co 

3 

at 

l 

1  im 

p- 

4. 

Vkt 

VI 

a 

2  2 

5  v 

<  “ 

vc  ~ 

. 

•u  1 

:  i 

!  X 

VI 

: 

> 

*• 

r 

a.1 

3  < 

,-  X 

r  x 

< 

as  2 

*  i 

"  ' 

c 

t 

1  J 

-J 

s 

=1 

a 

x  ►- 

P-* 

X 

x  at  < 

I  • 

J  ! 

-J 

X 

z 

J 

Zt 

^1 

-i  < 

■< 

<  < 

< 

VC 

i  x  a 

3  1 

«,  i 

< 

c 

< 

* 

o 

i 

«  < 

<  L. 

1  *  X  U  a 

-j  a 

X 

1  <  3 

<3 

z  ; 

1 

2  | 

2 

CL 

£ 

♦ 

O  u 

o  - 

X  c 

^4  C 

<  x 

> 

V 

|  V 

P 

1  | 

P- 

C3 

V* 

CD 

X’ 

—  a 

< 

5 

O  i 

< 

r 

wc 

z  ; 

i 

I 

VC  | 

v 

I 

25 

irt 

C 

•-J 

oc  a 

a  C 

X  MJ  "  M 

—  X 

XX*  vi 

c 

,  * 

« 

< 

a- 

< 

=L 

< 

X. 

1  c  c 

c  ► 

C  X  ^  * 

G  c! 

X 

i 

L  X 

*i 

G  i 

X 

w 

r— • 

u 

SL 

JJ  x. 

H*  ►“ 

k»  V 

:  </< 

w 

-0 

• 

l 

z  : 

1  : 

VC 

3 

>- 

VC 

1 

3  » 

v.  / 

V*.  •—  H  VI  •• 

vc 

4M 

-*•  2 

! 

«  ' 

f  i 

f 

1 

li 

2  i 

M  «<4 

•“  Z  v» 

X 

T*  ^ 

G 

H  1 

x  • 

jf  ■  * 

1 

» 

J 

X 

MJ 

—  X 

X  X 

X 

—  i 

1  1 

• 

04  I  W 

! 

i  I 

r  * 

►* 

G 

!  i’ 

-J 

< 

< 

<3 

«  X 

• 

1 

• 

X 

41  £ 

1 

: 

i  ! 

1  I 

3 

* 

v» 

— • 

•# 

Z  ¥ 

• 

•  < 

X 

*  t 

• 

1  X  V  v 

! 

1  1 

_i 

!c 

41 

K 

X. 

~  C/l 

4-*  f\ 

J 

5 

<  i- 

•w 

1  I  3 

1  - 

' 

x 

f  • 

I  [ 

X 

<* 

■ 

* 

i 

L‘  — 

i 

’ 

1  1 

r* 

1  1 

U^UUuO 


cuuuduowUud 


iJw 

\ 


D  35 


uUouuuuJ 


* 


I if. 


t  ; 


OC  O' 
O'  O' 
0>  O' 


c  c 
c  c 

r-si  eg 


li 


jS 

i*r 

D 


< 

A 


>*  c 
o 
c  c 

Cm 

c 

o 

c 

o 


o  c  e  c  c  c 

<■  i^r  £  c  ~ 

3  C'C  c  c  c 

CC'OOOC 
3  0.0  O  C  O 
3  DO  CCD 


>  X 


VI" 

(X  ••  Jx. 

*s,\ 

1  * 
••  I  *• 

•I-  s|n 

^  «  U  ^x 


< 
p-H  ’ 
<\ 
H  * 

X  . 


—  ,z 
»  —» 

*  - 
w  ^  •» 

X  «  3C 

►*  :i  r 

—  *  • 
23-. 
XX* 


I—  j 
°; 
s! 

yi 

&  ; 

IH  ‘ 

2: 

•  ~ 


*-p  *  i 

£ 

< 

O 


c  c  ;c  o  Sc  o 

|f\ 

c  !c  c  o  c 
c  c  c  c  c 
CCiCC’CC 

pc.cooc 

C  C  D  CO  C 


v>  t> 


yi  >s 
— '  !  • 
<  UJ 
IL  ca 
•  J 
►  N  < 
Ai  U. 

3  ‘* 
-J 

.  2.  ■V 
*  ( 
-J 

•*  C 


u. 

Bo- 

■-  y  »  x 


i 


Is  I  >.  u. 

Cj-3-‘ 


Ai  ^ 

2  ;v* 

if 


—  H  o 

w»  ^*  <  > 

—  *T  O 

-J  X  —  < 


ucu  — 


mJ  Ll 


t/1  * 

-i  * 

<  c 

x  :c 

•  c 

n  *! 

• 

C 
►-  V 

Z  x* 
•»  < 
i  x 
V 

<  < 

>-  X 

<  4 

D  C 


coo 

J  O  -H 

-  A.  OC 

\z  c  c  c 

C  C  C  i 

occ 

o  o  c  .> 


A.  cr, 
A,  A 

c  c 


pep 
v  < 
ev  <\ 
o  !o 


E  Dio. 


!* 

ur 


a: 

vo 


2 

uT 

-j 
< 
o 

VO 

>  • 

w  a.  V*.  -  ;T 

-<  ;  u  ..■  v 
*  •*  x  a 
c  »-  U  » 

.  •  •  *-  VO  1 

X  A  A  V.  -w 

—  rs  r  « 

*  ;c  -  a 

*“  »j  •  *  A  Z 

z  a  .%  z  »' 

u;  jr  -4  r  * 

I 


>  -« 

.  £  7 

**  c 

►  .1  ^ 
X  / 
x  ;  • 

A  ... 

>  c.: 


=  feg 

(V  '»%  f'l 

c  b  c 
c  'o  c 


,< 

c 


-  g 

G 

PI 


<7 

9  m~ 

is*  o 


Vi  a 

"  v. 

V  < 
..  x 

»"•  I 

W  A  . 
<  ' 


z  x  • 

X  K 


cope 
rv  k  iv 
m.  cr.  'x<>  cr 

c  c  jc  c 
e  3  c  c 


:occt 


D  C  C  D  P  D 


pee 
>  <r 

fp.  jpn 

pep 
^  c  c 
c  c 
p  o  p  c  c 
p 


P  ! 


X-  *v 

:  a 


a  . ; 

»•  •  i 


I 


x 


y.  ] 

A 


-  U-  V  ‘ 

A.  t'  •• 

•  <  V 
C  u  *» 

w  • 

A  j' 

i>4?< 

,  ►"  *  X* 

•  —  l  ►- 

,  i  ;  l~ 

i:  >  x  * 


41 


COO 

o  c 

c*  ^ 
c  c  c 
c  c  c 

CCD 
CCD 
O  O  O 


D  «l  v 
.>  >  > 
-II 

*-  A  A 

i  x 
D  >  > 
u 


£  -  > 
vt  V  v. 

;r*  r  t 

o  *•  vk 

.*  D  > 
D  a 

A  x  A 

;:  d  > 


O  I 
po  , 


h3; 

p  D  • 

5  A 

-li- 

V  X  i 
>  •?  -< 
D  "  • 

w  > 
A  A  Z 
l  C  I 
>  r  — 


\ 


4 


r 


i  ,i 


O  cl 

«T 
3>  e* 
o  c 
c  c 
s  c 
|o  c 

3  O 


O  CO  c 
If  7~'3'  O' 

c  c|e  c 
c  cte  s 
c  ole  © 

c  CO  © 

©  o  ©  © 


C  Cto  etc 
,c  0:0  o'c  c 

O  ©;©  ©I©  O 

to  ©;©  ©o  c 


jc  o  c  o'c  c 
c  ■**?N  t*  <  A 
C  c:c  CIO  c 


» 


« 


tt  A» 


» 

<MJJ 
-  © 
r 


cl> 

il 


C  CIO  c  c  c 
<  X  0s  c  ^ 


c  c  :c  c  c  c 


C  O'C  c  c 
SCjCCCC 


a  ,u. 
\> 


2  Z 


SIS 


li 


.  I 

%  r. 


a. 

12  1 

•  X 
*-  '< 
•  '  ,* 
X  > 

•  A 
w  K 
x  t 


X  x 

>  X 
I  I 
X  X 
<  < 
A  X 

►  X 


U. 
V 
X  £ 
>  — 

*.5 


AAA 


C  O 

Oi  r.  |sr  vr  r- 


C  C  ic  c  c  c 

bocccc 

COjCOCC 

c  o  b  c  ;c  o 


E 

z 


X 

X 


A  A 
£  > 


O  W 
\  / 
A  I 
2 


X  > 

>  z 

C  X 

*-  D  £  * 

4  •  A  ; 

C  A  —  ~ 


-A  A 

u 


bcco'cc 

•X  C  C  A  tr\ 
re  a 

'^4  '*»<  ^4 

bc'cccc 
b  o  c  c  c  c 


be 


_ _ bccccc 

kr  IT  ;*  P-  ix  S'  C  -«  A.  I* 

H  A'  a  A  ;a  A  tr\  tr\  rr 


Ecbc'cc 


p  o  s  o  c  o 


k  I  - 


lr 

>  r* 

tfl  CM 
*  ' 

» c 

N  >- 


—  A 

il-1 
;  t 
r  \ 
*  ,n 


►  c 

•  c  W 


\ 


I 


A*  A 
A  V 

15 

I  t 

-  .  -  XX 

A  ~  A  D  A  ~  — 

1  i  —  V*  —  x  X  £ 

£  1  £  K  ► 

w  -  A  A  ► 


> 

A 


Iff-  cfe 


c 

c  c  c  c 

C  C  P  o 

coco 


I 


beccico 

!§§§§:§§ 
o  c  c  :c  c 


(A 

Z  A 


►  IA 

•  t 

x  C 

*  h 

r 


U  Q 
w  C 


a 

IA  ** 


X  A 

X  c 


<fi  •- 


[CO 

,.*A  rr. 


fccoccc 

KC  A-  CC  ^  C  — 
tfT  <A  fT,  tn  < 
t— •  «— 

bccccc 


pccocccc. 

bOCC'CCCC 


I 


pococopo 


K 

o 


u-  —  a 

^  i  •• 


It. 


<  ♦ 

U»  4  ►“ 


43 


p  o ; 
rC  tf\ 
P*  * 


o 

A  ►- 

ro  i 
c 
c  c 


b 


Ss 

c 

~  Z 


3  U~ 


•  c 


li 


jK. 

•  sc 

Ai  c 
X  c 


«  r*  i  £ 

j:  • 

zz;  ^ 

-*  —  I  > 

ts  2  •  O 
:  x  >*  c  * 

*  I  I  •  • 

*  —  ►'  w 

- -  c 


X 

2 

V 

X 

> 


►  X 
X  w 


< 
••  X 


H‘'»A  K 

—  .*t  .^1  < 


—  •  S  3 


o 

Al 


c 

A* 

M 


<  a  a  : :  — 

—  < 

<  T  —  *  »~ 

-  ■ —  •-  x  ■ —  _  -  ■<  -  i 

0  a  —  a. 

£  £  c  ;  *  z 

C  0  w  X  4  A 

0  u.  I  ecu  J.  3.  t  1C  *  C 

i 

-  ~  w  r.  -  - 

C  —  y  x  >■  X 

w  —  s  x  o  —  -  v;  *  ■*  it  : 

t 

X  D 

c 

i  1  '  *  1 

3  —  A  f  I. 

A* 

A» 

IT.  ^ 

5  c  ^  ,  * 

K 

■  « 

N  IV 

A  A 

.  , 

J  t*.  t* :  —  *\  .  •-  ^ 

•  ^ 

:  i  ■  ■  ! 

I  — 

2  D  -  u  I 

i.  N  A  A  2  » 

r  r  —  ^  5  o 

T  't  w  w- 

i  u 


-  N  X 

s  r  • 
:£> 

z  -f  X 


M  c 


r  — 


A  A  — 
X  >  A 
O  X  "1 


A 


^rvAvkSm 


Option  NMP  Listing 


(Main) 

Subroutines 


READ 

STAT 

TTEST 

STAND 

CORREL 

OUT 

ARRAY 

MINV 

GMPRD 

GMSUB 

TRIANG 

FLOW 

CLT24 


XYSPA 


'h 

■1 


SimKUUTINF  TRIANG--  ENTER  C  HATH  IX,  SOLVE  FOR  TRIANGULAR  HATRI 


71ho4hSITI 


pm  m 

o  o 

PM  PM 


3 


vp 

c  c 

PM  C\. 


—  V> 

z 

<’  —  ,  • 
**.  o  — 

«C  •  N 

C;-C  *x 

U.  *- 
~  IN.  — 
33  -*  A.  V* 

iP' >-  r 
~  -  3 

C  JC  ^  - 
*  r  r 
m  :  c  c 
x  x  g?  - 


<*  * 


pi  m;*  ^  ^  r- 


< 

,3 

pi  at  * 

>  ac 

-  .  *  < 

—  3  •  x 

‘-  'r.  -  > 

Vi,—  H  — 

—  K  X  — 

— ;c  <  z 
^■x  *■  * 
•  —  >  c 

■x  «  3  — 

p*  2,  —  — 

—  X  X 

a  z  Z  £ 

X  x  z 


u 

?' 


2 

< 

i  *" 

j  VI 

< 

X 

X  Vi  V 
k-  <  — 

—  X  X-  • 
vi  >  v  jf 
Z  ZZ  < 

»  •*  •  x 

—  —  —  > 


»•!  *•  <  «  X 

r  x'  -  -  x  - 

—  Z  /  y  wi  .. 


X  (X  C  *H  PJ  fp 
«■*  —  ;PM  Pi  Pi  PJ 


hf  V>  >C  h*  X  a 
KM  PM  Pi  Pw  Pi  P. 


Pw  Pi 

Pi  Pg 

•N  PW 

PM  Pw 

Pi  p. 

'%tv 

*Si  Pw 

Pw  p, 

pw  rw 

%  <y< 

• 

c 

X 

■ 

WO 

X 

2 

» 

< 

0£ 

't 

1  | 

O 

• 

1 

• 

PW 

i 

X 

•“ 

1 

i 

X 

■ 

z 

IP 

— 

< 

Pi 

i  j 

l A 

oc 

■  f 

< 

o 

o 

X 

1  1 

Ui 

• 

» 

< 

00 

o 

< 

> 

z 

►- 

X 

< 

* 

♦ 

c 

II 

• 

v» 

h* 

3 

» 

h~ 

< 

Vi 

U.- 

*r 

Ui 

< 

X 

PM 

3 

V) 

UJ 

VJ 

M 

y* 

«k 

c 

z 

•t 

UJ 

ac 

v: 

+—  j 

A 

3 

1  * 

< 

Z 

\A 

> 

< 

-j  ! 

< 

< 

< 

X 

<  i 

U) 

> 

Ui 

o 

U  | 

v> 

> 

z 

+*  i 

•X 

> 

, . 

< 

* 

» 

H- 

Ui 

< 

D  ! 

at 

— 

Hi 

X 

- 

>-  ! 

< 

-J 

X 

g> 

wo  j 

Ui 

< 

* 

PM 

|> 

Irt 

WO 

X  | 

c 

i 

— 

< 

X 

■  i 

p  —  .Pi  m  vr 

rn  p"i  ro  rP  m  m 

K\  pm  [fM  Pw  pm  pm 


Xii 

vi ! 


wO  : 

Z  ; 

*  i 

>  1 
Z  : 


U‘ 

Z 


u. 

C 


z 

Ui 

X 

o 

X 


kC  P  ®  O'  c  ^ 
p  m  <«i  pi  !g-  wi¬ 
ry  <\  <M  PiiPi  PM 


c 
•  c  c 
Oj  •  • 
HOC 

M 


|N  m  :>*  tP  >C  f*» 
<  ^  [-tf  <■  wj* 

[Pi  pm  j1 fM  pm  ,pj  pm 


r 

(V 

■< 


WP  5 


t  ! 


s 


X  O'  o 
•T  -t  ^ 

i<M  <M  ,<v 


*  r* 

nr. 


C 

JJ 

b 

Vi 

X 

a. 

VI 

< 

J 

C 

■  ^ 

z 

it 

tl 

« 

\A 

• 

• 

< 

'X 

xK. 

Vt{  • 

u 

u 

< 

M 

»■ 

x> 

JZ 

Vi 

vs 

- 

vx 

> 

•  * 

<  V 

« 

X 

v> 

X 

•Vi 

\xi 

< 

< 

wn 

V 

> 

— 

-  i.»i 

Ui  < 

V 

x 

o 

2 

& 

X 

Z 

wO 

'JJ 

X' 

z 

r. 

Z 

«k 

c  {  • 

Vi  X! 

x 

1  * 

* 

• 

— 

V* 

VI 

m 

• 

X 

“i,N 

«-|Vi 

• 

o 

L 

ML' 

«•* 

• 

"• 

•“* 

p- 

r  s 

X 

X 

z 

H  -IK- 

n 

*4. 

V' 

U  i» 

It 

M 

— 

• 

X  1  • 

►- 

< 

u  • 

V 

u. 

u 

•  X 

V 

X 

•* 

•r  y 

M 

C  *-  x 

•• 

«— 

✓ 

< 

— 

k- 

►» 

—  ■  *» 

< 

< 

•  V 

M 

•  h- 

s 

Pi 

yt  * 

<  x  x 
;x  x  x 
\A  ♦  ♦ 

X  X 

hx  x  x 


Iwo  h-  *~ 
\A  J\ 

I*  -  - 

< - 

>  X  3 
1*  «  £ 
H  z  Z 
—  <  < 
xx  x  x 


W  —  X  —  H 

- - -  -  v'  /  V 


—  *  'iJ  -  i 


5  v»  —  • 

►.  m 


r  r  :  i  < 


c  ‘ 
-  ?  g  u  y 


Z 

vt 

IT 

_ 

X 

w 

3  •»* 

— 

X' 

— 

»• 

— 

w. 

* 

X 

1 

X 

Z 

X 

z 

z 

P* 

X 

V. 

z 

< 

3 

X 

2 

X 

Vi 

X 

u 

V 

•M 

w 

i 

j 

p. 

N 

\ 

Pw 

<  *- 

•» 

X 

-f 

X 

* 

X 

^  c 

X 

z 

< 

M 

< 

X 

< 

U 

M 

• 

M 

X 

< 

< 

< 

X  z 

< 

X 

X 

4-  C 

i 

> 

“ 

C 

X 

“ 

•.  X 

X 

* 

X 

X 

> 

:t 

X  x 

V 

X 

‘X 

V 

- 

“ 

C 

w  r 

X 

V 

- 

- 

- 

X 

C 

X  J 

X 

X 

K 

X 

* 

“ 

U  vJ 

i 

I 


oc  u 


MJ  P' 


r/r.RANsor  isttetk 


1. 


cm  c«i|  gr  v* 
o  c;c  c 

PO  fT,:  (‘T  fC| 


O-JO 

ft-1  H 


Vl>  ¥ 

X  ¥ 


Vi 

2 
l  * 


>-  i*  .i 
X '  •—  ft  x 
^ .  V  ft-  ft- 

—  ft-  *  — 

V*  —  ¥  Vi 

—  ¥  ¥  — 
* 

*  X  *3  T 
< 


'C  Is-.  rc  0-'  O  «-• 
C  C  C  C  -* 
cci  cr:  m  r«vro  <rt 


|(M  nv*^  ^ :*c  o- 

<— •  »— I  —I  «— 4  «— < 

on  (tyro  (T\  (r,  *<y 


X  *  i  t  ’  w  ~ 
r  ^  2  .  «•*  »-  *— 


’ 

, 

«r 

ft 

—  «Nl 

a: 

<  - 

>  Vi 

X  CM 

ft  ft- 

Z  ¥ 

<  < 

UJ  > 

I  X 

X  • 

ft  —ft 

» 

<  >f 

i 

ft-  * 

<  CM 

¥ 

2  —4 

VI 

••  • 

UJ  — 

1 

X  Vi 

SA  ¥i 

ft-  CM 

>v  ¥' 

l 

— 

vi  Z 

¥  Xt 

z  < 

¥  ►-: 

, 

•  UJ 

ft  •"  ; 

X 

3t  X 

UJ  w 

3 

>  X 

ft-  —l 

•— 

X  • 

ft*  ftft 

H* 

» 

V5  X 

2 

1 

V  CM 

<  - 

V 

*i  • 

X 

¥ 

at 

Vi 

vi  O 

c 

V» 

»— 

a. 

Z  O 

VI 

z  + 

Vi 

V. 

•»  — < 

¥ 

<  — ' 

—• 

•ft 

••  » 

X 

‘at  x 

¥ 

—  *s 

Vi 

•J  * 

¥ 

—  PM 

e  _ 

X 

••• 

2 

V 

a 

kJ 

2  < 

X 

►- 

¥  Vi 

¥ 

Vi 

<  ft- 

¥  — 

4J 

•ft 

ft-  < 

o 

¥» 

«  22 

^  Z 

V?  2 

X 

x  ^ 

’X 

*  ft- 

•ft 

X 

►-  m 

—  u 

x  X 

ft- 

Mr 

ft-  ^ 

a. 

II  *— 

X  c 

— • 

< 

Vi  ¥ 

¥  ¥ 

—  — 

2 

'X! 

•-  ¥ 

*— 

►- 

¥  ¥ 

—  Vi 

¥ 

ft-  w 

— 

«• 

r 

:  z 

< 

X 

X  — 

/ 

2 

c  — 

X 

c  *> 

2 

ft-  V 

X 

—  V 

z, 

*  r 

X 

’X 

ft*  ft- 

it 

«4 

•ft 

■ft*  * 

^  — • 

< 

Vi  X 

✓ 

X 

x  ■* 

r 

*  2 

ft* 

w 

»—■  *— 

X 

■  ^ 

¥  a. 

A 

V 

CC  0s  ;c  —  og  <*i 

^  — /cm  <\  rvj  rs, 

fTj  (<\  fO  fO 


I*  ir  e  f-  x  <r 
KM  <\i  CM  CM  CM*  CM 
Wi  m  .  <%  <r.  r**i  rr, 


C  — 5V«J  V*  i£ 

,*0  (T. ,  ,"0  «r. '  re  Lv. 


O' 

cr 

<M 

O 

m 

• 

fVJ 

—  I 

vr 

< 

X 

V> 


5 


Vi 

¥  • 

•  'JC  X 
Vi  <>  Vi 

2Z2 


i  21 

:  <! 


Cl 

-It 
H  , 

v»i 

2; 

Vl| 


rr,  m  re 


y ' 

< 

X 

-  ft  rr  >• 

—  — ■*  tet  ft  te- 

*  I H  N  CM  * 
x  i  v*  •  a 

—  »x  JJ  —  — 

V  *>  w*  .%  v-  «• 

—  —  —  l  —  ^ 

— *  w- 

:  I  *  i5  i 


X  *¥  X 
■212:: 


x  ,l  y 

—  —  > 

•  W*  V  X 
••  ‘r-  V. 

<<  — 
w.  v»'  ft. 
<  <  - 


w 

2 

< 

►«• 

Z 

w 

X 

3 ' 

X 

< 

< 


-*  r  c 

,  m  rr  !>f 


V. 

< 

X 

Vi 

•ft 

• 

-— 

lu- 

ft- 

vi 

V) 

Vi 

•* 

c 

z 

te 

. 

— » 

< 

*ft 

X 

•- 

IT 

Vi 

X 

X  CM 

< 

X 

• 

•x 

► 

X 

1 

ft-  IT 

Vi 

>  CM 

te¬ 

Vi 

2 

— 

¬ 

< 

o  a: 

ar 

x 

* 

X 

< 

v: 

<  ft- 

X 

— 

a. 

> 

> 

• 

ifl  c 

•- 

¥ 

> 

» 

< 

X 

c 

u 

—• 

ft- 

> 

•  Vi 

H- 

te- 

<  CM 

Vi 

- 

*- 

- 

ft- • 

b 

<  Vi 

— 

ft- 

j.  ’ 

3  ru 

< 

< : 

• 

— - 

►- 

Vi 

>1 

X 

¥ 

1 

< 

X 

X 

1 

2 

-— 

oi 

te-  v> 

1 

w 

< 

X' 

Vi  > 

■  1 

V— 

ft- 

2  O 

1 

¥ 

< 

•ft 

- 

j  i 

—  1 

O 

3  — 

<2 

X 

«  i 

t*, 

v. 

— -  te¬ 

5; 

*- 

ft 

M 

ll  te. 

fV 

X1 

•• 

-ft-  ft 

2 ' 

< 

-  1 

2  C 

- 

VI 

Vi  Vi 

xi 

X 

-  . 

•  • 

2  ! 

¥ 

<  < 

<f 

<£  C  1 

3  O 

< 

X 

JJ  X 

•* 

2  0! 

II  u 

V 

vi  v. 

Vi! 

*• 

-  i 

—  ] 

Z 

ft- 

—  — 

i 

» 

JJ  X 

X  X 

O 

XI 

•¥.  '«£  ft- 
n  'f  'f  «:  -o  >j- 


X  a  c 

<■  Vi 

r*i  ro 


n. 

X 


ouu  — 


u  <  v  ~  b  x 

C  <  H  W  w  X  . 
a  x  a.  u  —  — 

*  3  v  ►-  *•  ^  vi 

ft-  o  v*  w'  —  — . 

2w*  r  —  ■>  >;  2 

xx  a.  * 

XX  M  —  — 

£  *  >  > 

.  v  x  x  z  r  j  u 


X 

X 

x 


x 


wn !  v) 

XXX 


H|K  » 
X  lX  ¥. 
*-  <  < 

—  'Xi  X‘ 
vt  >.  V 


\  M  \ 

r  x  c 


I  —  U  U 


J 


?  ;'5'2 


•-«  cm  ml*  ir>  < 
u%  ir>  irjir  *.  tr 

«rv  fO  (n|*o 


f  i  !  i  i  ; 

r-«  T'C  «■*  .*vj  ro,^  0s 

x^tT’C  x:  -c  x  >c  >c  x  ^  <c  >cp**fs.l^rfcp^p*-p*-^.p»-fs» 
**■ j ^  c*i  cn  m  rn  ,*o  .•f'  (r,  (Ti  f\  rn  «r>  (<\ 


—•  fsi  c*\  i**  x*  -C  p-*  cl  o*  c  **  m  It  tr  Lc  r-  lx  o*  cr 

X  X  X  cc  cc  •::  X  CC  X  O  O'  T  O'  O'  O'  'O'  O'  O'  O'  c 

.-r.  pn  ,r«  ,*r,  jrv  <  ^  p  m  fn  rr,  ^  m  r*  m  m  m/o  m  -t 

II  I  II 


;  vj  •"  —  » 

i  «*.  iA  X*  * 

b-  JT  —  , 

‘X  *  •  * 


“  b*  *-  . 

.  —  *-.  Jj 

I  w  *  - 
>  <  ^  X 

1  *  — 
>  *-  r 
z  <  z 

U  £  *J 


U  O  ^  f>. « 


b>  f  *-  «-*/»•. 

• •  *>*  ^  ^  U.’ 

•>.*"*  *  x  :t 

u.  —  Z  C  w  X  •-* 

O  5  ►». 

>  -  >  •  «r 

r  z:  <  z  z  c  *  c 

^  w  r  l*  ?  'j  * » 4 


x  < 

u.  X 

K  ^ 
>  > ! 
O  O: 
u  o 

t)  II  I 
X  —  —  ! 
>•  u-  x  i 

~ , 
a- 

•  "in 


■  i  -  — 

x  <  r 

«  ‘i  a.  a 
—  »-  Z 
>  >  *- 


•n  -  c  u> 

‘  O,  K 


►*  C  I 

«2  k 

«  •=  !• 


\ 

Jl 


turn *«Whb-  ? 


IT 

XlX**-  X 
1/vJwA  i^'u^  if\ 


Z 

C: 

Ol 

UJ 

*! 


*1 

v»{ 


O! 

© 

•i1 


► 

x, 


*c 

? 

•** 

w 

4 

ar 

i 

e 


£ 

'A  \f\.t 

/>  (/V  I  X 


O'  C 

jT’C  •£ 

*A'A  iT 


u»a  <\ 
—1< 
✓>|a. 

5!§ 

^iU 
'  * 

WT  — * 
— '  iTi 

C;AJ 
►-  - 

•►ixi 

*>% 

••x 

<v(o. 


X 


X  J- 


•  > 

»rjc 

u  * 
**  - 
v/I.X' 

•A 

*;* 

* 

NiZ 
|—  < 
X 

H  *  C 
*  * 
VI-'*' 
<  »  ^ 
,*J  <•  * 

is*  > 

Z  *  <\ 

«  —  * 
>  +  a 
Z  *X 
*  A.  < 

u  -m  at 


>■}—  ^  * 
*tV»  N  <•» 

iT  ?  *  ^ 
*  X  • 

*  X  «  a 

r  x  j  -  in 

C'*»  i  - 
»  i 

U  1  M 

*  JJX 
"•  •  C  H 

w  3  ’«£ 

*  x  *  • 
»>•*■** 
£  *  s  * 

¥**  '*•  r.  •* 


-  f 


W 


i  < 

|*C  <  -C  <C  ^  X 
*  «T  <A  *T  u*n  T- 


4 


Z 

c 


0*  ©  —  M  rr. 
•£  A-  fw'A-  >. 
S<  if\  A.  X 


'*  •/*  *c  r-  x* 
is*  p*  a- 

1^  iT  J1'  J*  X 


*  £: 
I  >• 

!  ij 

.  * 


c  i 

X  ' 

Si 

wu  j 

*i 


I 

I 

5 

\  < 

4*  r 


z\  r 

—  5 

"  **,  * 
»*  St. :  ~ 


+  <  M- 

—  *  3 

—  *  a  .j 
*  i  *  i 

a-*  ^ 


< 

O 


X  . 
X  * 
o 

N 
•  • 

X  ■• 

o 

N 

3 

I 


•  >f  :  / 

\  • -<  ^  *( 

>  ^  ~  ! 
~  U.  *t  C  U  ‘ 

K  >«  £  ►  ft*  1 

--  a  —  - 

a.  t  :  ;  .  ■ 


X 


PA 

v> 

r* 

v*. 

> 

Lj 

'< 

* 

.< 

'VZ 
N.  — 
wn  i 
-  C 

i/!  u. 

>• 

— j  z: 

< 

Z  UJ 
<  u 
< 
c  « 

X  w- 

o  c 

•xi  UJ 

*« 


>>  > 


^  CC 
u  Uj 
Z 
<  u. 

3-  O 
u-  , 

*  U. 

u;  o  . 

O  ' 

X  trt 

*  X 
A»  41 

•  W 

K  >  «  < 
O  4^  WJ 

a  ©  r  a 
O  o  Z  ,s 


—  ✓  V* 

r  v  <  -« 


u.  ^  v* 


& 
s 

i  s 

~  V-  ? 

«  r  • 

•  -  f 


<\i  r\ 

X  X 

X  pft  X  O' 

e« 

a  m 

^  .r1 

X  X 

X  x 

X  X  X  x  :o  a* 

>  O' 

O'  O' 

x*\  x- 

X  X 

ft*.  X  IT  x 

1 

1 

x 

X 

A  x 

1 

c 

N 

1  • 

* 

> 

UJ 

E 

i 

..i 

, 

wft 

iK 

X 

N 

p  , 

X 

< 

1 

pc 

UJ 

UJ 

X 

•— 

Wo 

•» 

X 

UJ 

fT 

►- 

r* 

w 

X 

rg 

ft- 

r* 

"*•  — ' 

r* 

^  X 

1 

I* 

Ift  < 

c 

>  1 

I* 

<  x ; 

l/> 

jj  • ! 

X 

A  -ft  ! 

'r 

-  * 

1 

1  ft 

ft  - 

X 

u  (A 

I 

N 

-  « 

ftft 

-  UJ 

(ft 

A  lA 

i 

—  —  1 

« 

—  -» j 

N 

UJ 

3c  uj  (  1 

-> 

-  »;  I 

a  h* 

*■ 

4j  «**  ! 

Ui 

C  X  j 

X  *  i 

“*  1 

X.  *  1 

*• 

H  x  i 

^  w  ,  * 

a, 

«  rw  « 

►ft  X  i 

-»  —  ‘ 

/  < 

o 

ft  XJ 

ac 

ift* 

£• 

s; 

c  - 


UJ 


^  x*  A, 
►  f  S^ 


«C  -  — 

c  ^  c 


NMB  it 


O'  >  S'  o 
1  it'  'C 

I 

I 


-H  5 
•  > 

^  Ul  X-'  U 

Zmrz  Ji 

~  w  ^  * 

«  H  —  * 

i,  k*  •  / 

•  V*  A  X 


o 
•  • 

uu 


•?  r. 

A  w 


4  w  < 


X  * 

—  3 


X  *- 

;  x 


x  ; 


/  S 

—  a;  C  —  X 

^  >  4.  X-  • 

-  y  f-  -  - 

s  -  •  r 

T  %  4  /i  W 

z  u  ?-  —  u 


e* 

V 


c  \ 

X  • 

—  J- 


►c  > 

■  -«  x 


~  X 
w,  3  • 


w>  w 


*-» 


^  ^  [  * 


I 


<  ik  i«r (  i « > *ST3TT*»T 

I  IIIKM 


OPTION  B  LISTING 


MAIN 
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Overlay  2  Routines 

CHAIN6  .  (Main) 

Overlay  3  Routines 
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EIMTHATHJNS  of  RESERVO  IKS  with  nonzero  starting 
centrations  are  used  for  EACH  RON. 
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SCALES  FLOWS  AT  GAGING  STATIONS  TO  INFLOWS  INTO  EACH  COOKO 
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